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THE  USE  OF  TWO-STATION  lONOSONDE  DATA  FOR  lONCAP/ASAPS 
MODEL  VALffiATION  ON  TRANSPOLAR  TWO-HOP, 
HIGH-FREQUENCY  SKYWAVE  RADIO  LINKS 

INTRODUCTION 

BACKGROUND 

An  investigation*  of  one-hop,  high-frequency  (HF)  trans- Arctic  skywave  propagation  was 
conducted  through  a  comparison  of  Scott  Base  (Antarctic)  ionosonde  measured  data  with  a  stand¬ 
alone  prediction  system  (ASAPS)  from  the  Australian  Ionospheric  Prediction  Service  and  with 
lONCAP  from  the  Institute  for  Telecommunication  Services  of  the  U.S.  Department  of  Commerce. 
During  this  smdy,  a  FORTRAN  program  called  lONOLINK  (using  the  modem  format  for 
ionosonde  data)  was  produced  for  the  creation  of  oblique  frequency  availability. 

The  results  of  the  investigation  showed  that  both  the  lONCAP  and  ASAPS  models  predicted 
the  single-hq>  pn^agation-dependoit  frequency  of  cq)timum  transmission  (FOT),  maximum 
usable  frequency  (MUF),  and  highest  propaption  frequency  (HPF) — ^frrequencies  not  exceeded 
10, 50,  and  90  percent  of  the  time,  respectively — ^within  2  MHz  for  a  wide  range  of  sunspot 
numbers  and  solar  zenith  angles.  (No  attempt  was  made  to  determine  the  lowest  usable  frequency 
(LUF)  values  because  th^  are  not  only  dependent  mi  propagation  but  also  on  the  radio-link  power 
budget)  Hie  close  agreement  of  model  predictions  with  the  measurements  can  reasonably  be 
explained  by  the  probable  crnitributions  of  the  Scott  Base  ionosonde  data  to  the  development  of 
itmospheric  electron  density  profiles  used  in  both  the  ASAPS  and  lONCAP  programs. 

SCOPE 

In  this  report  the  one-hop  results  frrom  the  earlier  study  have  been  extended  to  the  two-hop 
case  using  ionosmide  data  simultaneously  acquired  frnom  two  high-latitude  (Arctic)  stations. 
Altiiough  ionosonde  data  frnom  a  single  station  permits  usable  oblique  firequency  extrapolation  fixim 
tile  vertical  incidence  measurements  to  any  single-hq>  range,  the  use  of  two-station  ionosonde  data 
restricts  the  results  to  a  single-hop  range  for  each  pair  of  imiosondes  used.  In  addition,  the 
ionosonde  locations  must  be  chosen  so  that  their  separatimi  distance  excludes  one-hop 


*  Submarine  Electromagnetic  Systems  Dqiaiiment,  High  Frequency  (HF)  and  Meteor  Burst  Communications  in 
a  Polar  Environment,  NUWC-NPT  Technical  Document  10,375,  Naval  Undersea  Warfare  Center  Detachment, 
New  Umdon,  CT.  30  September  1993. 
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propagation,  while  the  dominant  mode  is  limited  to  two-hop  propagation  (i.e.,  the  range  should 
exceed  about  4000  km  but  not  be  much  greater  than  8000  km).  The  comparison  of  lONCAP  and 
ASAPS  predictions  with  these  two-h(^  frequency-availability  values  provides  a  mcae  accurate 
validation  of  the  model  for  complex  links  than  could  the  one-hop  cases  of  the  previous  study. 
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APPROACH 


lONOSONDE  DATA  ANALYSIS 
Northern  Latitude  lonosonde  Data 

lonosonde  data  measured  at  several  Russian  and  Canadian  ionosonde  stations  from  1991  were 
available  for  the  characterization  of  local  ionospheric  skywave  propagation  conditions.  The  data 
ctHisisted  of  monthly  tables  of  hourly  values  for  several  scaled  parameters  measured  from 
ionosonde  output  traces.  One  data  file  was  provided  for  each  scaled  parameter  for  one  or  mcno 
months  during  1991.  These  parameter  files  included  the  critical,  or  maximum,  frequency /v 
(MUF)  fcs*  which  the  ionosphere  would  reflect  a  vertically  launched  radiowave  firom  one  or  more 
icMiospheric  layers.  Ideally,  the  virtual  reflection  height  h’  would  also  be  provided  or,  at  a 
minimum,  the  scale  facuv  for  determining  the  basic  MUF  value  fca*  single-hq)  propagation  to  a 
fixed  distance.  The  F2,  FI,  and  £  layers  define  the  most  important  itmospheric  altitude  bands  used 
for  HF  skywave  radio  and  therefore  are  of  significant  importance  in  the  analysis  of  ionosonde 
parameters  as  depicted  in  figures  1  and  2.  For  the  F2  layer,  a  MUF  factor,  M3000F2  (Mp2\  is 
provided  to  convert  the  F2-layer  critical  frequency,  foF2  (where  o  defines  the  ordinary  wave),  into 
the  conesponding  MUF  values  for  a  single-hop  range  of  3(X)0  km.  Note  that  MUF  refers  to  the 
maximum  usaWe  oblique  fiequency  at  a  given  time,  not  the  50-percent  value  fix)m  the  distribution. 
At  this  point,  we  will  refer  to  this  "basic  MUF"  value  as  the  maximum  propagating  finequency 
extrapolated  to  a  specified  range  fix>m  the  midpath  ionosonde  data  and  let  the  "unmodified  MUF' 
refo*  to  the  SO-percentile  value. 

The  iixiosonde  data  for  each  month  and  scaled  parameter  were  presented  in  a  single  block  and 
smred  as  a  separate  file.  These  files  had  to  be  merged  to  create  a  single  file  containing  all  the 
parameters  scaled  for  each  of  the  available  months  in  order  to  produce  the  necessary  input  file. 
Selected  files  were  merged  to  create  input  files  for  the  IONOLNK2  computo-  program,  which  is  a 
FORTRAN  software  analysis  package  developed  to  determine  two-hop  finequency  availability 
using  scaled  ionosonde  data  fitxn  two  stations.  Icxiosmide  stations  woe  selected  to  provide  a  total 
two-hop  path  length  betweoi  4000  and  8(XX)  km  and  a  second  path  length  above  8(XX)  km  (the 
nominal  upper  bound  on  two-hc^  path  length).  In  each  file,  a  block  header  specified  the 
ionosonde  site,  the  date,  and  the  identity  of  the  scaled  paramet^  whose  values  immediately 
followed  the  header.  Each  scaled  parameter  in  a  block  consisted  of  five  character  positions — three 
positicms  carried  an  integer  number  and  the  other  two,  a  descriptive  and  a  qualifying  letter.  The 
scaled  fiequency  values  were  presented  in  the  data  file  in  units  of  \0f^,  MHz,  with  the  heights  given 
direcdy  in  kilometers  and  the  basic  MUF  factors  presented  as  \OOMf.  Many  hours  of  various 
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Figure  1,  The  ik-Value  Versus  Range  D 


Figure  2.  Two-Hop  Link  Geometry 
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months  and  years  showed  no  parameter  values  or  explanatoiy  letter.  A  complete  description  of  the 
data  format  is  provided  in  appendix  A. 

OBLIQUE  INCIDENCE 
One>Hop  Extrapolation 

Let  /v  be  the  ionostmde-reported  critical  frequency  for  a  vertically  incident  signal  fnxn  scnne 
ionospheric  layer  at  a  virtual  height  h\  The  corresponding  basic  MUF  value,)^,  for  single-hop 
oblique  propagation  over  a  link  with  great  circle  path  (GCP)  length  D  and  the  path  midpoint  located 
above  the  ionosonde  may  be  given  by 

fob=k(D)fySec(^o)  ,  (1) 


where  the  angle  of  incidence  of  the  signal  incittent  on  a  curved  ionospheric  layer  is 


0«-tan 


-1 


sin(gn/2) 


l+(/r7/?e)-cos(do/2)J 


(2) 


and 


Od 


(3) 


In  this  equation.  Re  is  the  earth's  radius  and  k(D)  is  the  secant  correction  factor  thc^t  approximates 
the  effects  of  ionospheric  curvature.  The  geometry  represented  in  these  expressions  is  depicted  in 
figure  3.  The  function  k(D)  is  plotted  in  figure  1. 

The  elevation  angle  0ei  is  measured  from  the  horizcm  to  the  virtual  reflection  point  and  provides 
a  means  of  estimating  the  usable  portion  of  an  anteiuia  pattern  during  link  design  studies.  From 
figure  3,  this  angle  can  be  determined  as 

e^=m-a-ej^/2 ,  (4) 


where 
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Figures.  The  HF  Link  Geometry 


a=:cos 


2rj.(RE  +  h') 


(5) 


and 


r;r  =  2Rs[l  -  cos(a^  /  2)](Re  +  h')  .  (6) 

When  the  virtual  height  is  not  available  in  the  icmostaide  data  base,  then  it  may  dther  be 
computed  from  the  corresponding  basic  MUF  factev,  Mf,  or  the  average  height  of  the  layer 
(e.g.,  325  km  for  the  F2  layer).  If  the  basic  MUF  factor  is  available,  then  the  virtual  height  may  be 
approximated  by 


h'=R^ 


“£ 


tan(^-) 


where  it  is  necessary  to  ctnnpute  tan(0o)  from 


(7) 
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These  equations  are  used  by  the  IONOLNK2  computer  program  (see  appendix  B)  to  determine  the 
critical  fiequencies  and  elevation  angles  fcH*  one-h(^  oblique  propagation  with  a  path  midpoint 
(reflection  point)  above  each  ionosonde  station  along  the  two-hop  path. 

For  each  hour  for  which  the  ionosonde  data  include  a  legitimate  scaled  value  (i.e.,  for  which  a 
value  is  provided),  the  IONOLNK2  program  first  computes  the  F2-layer  basic  MUF  value  for  each 
range  D  that  is  input  by  the  user  (i.e.,  fF2(^))-  Next,  it  computes  fpjiD)  and  fg^D).  Since  near- 
vertical  incidence  skywav*  (NVIS)  propagation  may  include  two-hop  propagation  fiom  the  E 
layer,  IONOLNK2  also  computes  the  two-hop  value  f2^D/2).  Finally,  the  sporadic  E  layer  (£,)  is 
considered  for  the  computation  of  the  one-hop  value/£,(D)  and  the  two-hop  value  f2EjiDI2).  The 
basic  MUF  value,  or  link  frequency /li  for  die  given  hour  is  then  determined  fiom 

fg  =  ^AAX.[fp2,fpi,f£,f2£,fEs>f2Es^' 

The  values  are  used  to  increment  cumulative  distribution  bins,  which  are  normalized  by  the  total 
number  of  usable  ionosonde  data  hours  in  the  chosen  sample  set. 

Two-Hop  Extrapolation 

Consider  the  two-hop  HF  propagation  geon^try  shown  in  figtue  2.  Two  ionosonde  stations 
are  shown  spaced  by  the  earth-centered  GCP  angle  Bs-  The  two  ionospheric  reflections  above 
these  stations  occur  at  virtual  reflection  heights  hi  and  /t2,  as  shown  in  the  figure.  The 
corresponding  "included"  angles  0i  and  respectively,  must  be  determined  to  compute  the 
desired  MUF  factors  (see  equation  (1)).  To  find  these  angles,  the  location  of  the  ground  reflection 
point  between  the  two  ionosonde  stations  must  first  be  found.  This  reflection  point  cannot  be 
determined  analytically,  but  may  be  solved  for  numerically  by  using  a  root-fiinding  algorithm. 
Root-finding  algorithms  begin  with  an  estimate  of  the  desired  root  (reflection  point)  and  then 
proceed  to  miiumize  some  oror  function  determined  firom  the  subsequent  root  estimates.  In  this 
case,  it  is  known  that  the  angle  of  incidence  {Or)  equals  the  angle  of  reflection.  Given  0\  and 
Bl-Og-  6i,  two  estimates  of  the  angle  of  ground  reflection  ( &ri  and  0:^2)  can  be  computed  and 
used  in  the  error  function  as  follows: 

«  =  «rl-«r2*  (10a) 
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whore 


with 

a-.  =  + /^  -  2r.R^  cos(0,.) 

and 


r,.  =  +  for  1  =  1,2. 


(10b) 


(10c) 


(lOd) 


The  root-finding  algorithm  attempts  to  select  6i  such  that  l£i  is  less  than  some  small  positive 
number.  Given  an  estimate  0\,  the  ground  distance  covered  by  the  total  two-hop  path  is  then  given 
by 


D^2(0i  +  02y*RE .  (11) 

Thus,  the  use  of  ionosonde  data  from  two  stations  to  characterize  the  overhead  conditions  for 
ionospheric  reflection  restricts  the  extrapolatitm  to  a  single  two-hop  length.  However,  the  actual 
endpoints  of  the  link  will  change  as  the  reflection  heights  change.  If  the  ionospheric  conditions  in 
die  vicinity  of  an  ionostmde  are  ctmsideied  identical  over  the  sky  volume  in  which  tl.e  correspond¬ 
ing  ionospheric  reflectitms  occur,  then  the  relative  frequency  availability  distribution  for  a  fixed- 
lengdi  link  should  be  insensitive  to  exact  geographic  location.  This  assumption  is  necessary  to 
justify  the  empirical  calculations  of  two-hop  frequency  availability  between  fixed  points  using  the 
lONCAP  or  ASAPS  computer  program. 

There  are  16  possible  two-hop  propagation  modes  considered  in  this  analysis:  F2-F2,  F2-F1, 
F2-E,  F2-Es,  F1-F2,  Fl-Fl,  Fl-E,  Fl-Es,  E-F2,  E-Fl,  E-E,  E-Es,  Es-F2,  Es-Fl,  Es-E,  and 
Es-Es,  where  die  hypten  (-)  signifies  ground  reflection.  The  IONOLNK2  program,  mentioned 
earlier,  was  designed  to  determine  the  nKxithly  two-hop  FOT,  MUF,  and  HPF  values  as  well  as 
frequency  and  mode  density  functions  over  tlM  entire  processed  time  interval.  This  program  inputs 
ionosonde  data  frtxn  two  stations,  determines  the  resulting  two-hop  oblique  propagation  geomeuy, 
and  then  calculates  the  basic  MUF  value  for  each  of  the  16  possible  modes.  Of  course,  links 
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exceeding  about  4000  km  will  not  experience  two-hop  E-layer  modes  because  of  the  lower 
reflection  heights  as  omipared  with  the  F-layer  modes.  IONOLNK2  then  determines  the  minimum 
of  the  resulting  basic  MUF  values  as  the  basic  two-hop  MUF  value.  This  value  is  used  to 
increment  die  appropriate  probability  bins.  The  10N0LNK2  computer  program  is  described  in 
appendix  A. 
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ANALYSIS  OF  RESULTS 


DIKSON  (RUSSIA)  AND  RESOLUTE  BAY  (CANADA) 

Scaled  ionosonde  parameter  files  from  the  Dikson  (No.  373, 72.5*  N  and  80.4*  E)  and  the 
Resolute  Bay  (No.  974, 74.7*  N  and  95.0*  W)  ionosonde  stations  were  used  to  produce  a 
7100-km  two-hop  HF  path.  Although  12  months  of  scaled  parameter  data  were  taken,  only  6 
months  of  Resolute  Bay  data  were  provided  for  this  study.  By  merging  the  appropriate  data  fries 
(see  appendix  A),  a  single  ionosonde  file  was  created  for  each  station  for  January  through  June 
1991.  Figure  4  is  a  plot  of  the  frequency  density  fuix^tion  covering  the  entire  6-month  period  for 
which  data  from  both  stations  were  available.  The  values  plotted  in  the  figure  represent  the 
percentage  of  days  during  the  6  months  when  simultaneous  dual-station  ionosonde  data  from 
Dikstm  and  Resolute  Bay  could  be  used  to  extrapolate  to  two-hop  oblique  propagation.  The  small 
number  of  usable  hours  is  due  to  lack  of  parameter  scaling  at  one  or  both  sites  and  to  the  dominant 
E-layer  reflections  that  prevent  adequate  two-hop  range.  The  month-t^-month  comparison  of 
icmosonde-deiived  frequency-availabiliQr  values  in  the  form  of  FOT,  MUF,  and  HPF  versus  time- 
of-day  are  presented  in  appendix  B  (pages  B-3  to  B-20).  Both  lONCAP  and  ASAPS  computer 
model  predictions  are  plotted  in  these  figures  along  with  the  ionosonde-extrapolated  results. 
Overall,  the  predicted  values  show  less  diurnal  variation  than  the  corresponding  measurement- 
derived  values,  as  shown  in  the  MUF  comparison  for  February  on  page  B-7. 


•160400 

>100-160 


UNMDtSAL  TWC  (HRS) 
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Hgure  4.  Six-Mmith  Frequency  Availability  for  Dikson  and  Resolute  Bay 
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CAPE  SCHMIDT  (RUSSIA)  AND  LOPARSKAYA  (RUSSIA) 

A  sec(xid  set  of  ionosonde  stations  was  selected  to  form  a  9100-km  two-hq>  trans- Arctic 
skywave  path  from  Cape  Schmidt  (No.  669, 68.8*  N  and  179.5*  E)  to  Loparskaya  (No.  168, 
68.0*  N  and  33.0*  E)  in  Russia.  A  full  12  months  of  ionosonde  data  from  both  stations  were 
available  for  empirical  extn^lation,  although  only  March,  July,  Septembo*,  and  Elecembo'  were 
chosen  for  die  month-by-month  axnparison  with  the  prediction.  The  longer  link  range 
significantly  reduced  the  number  of  hours  in  which  die  extrapolation  yielded  two-hop  propagation, 
as  shown  in  figure  5.  The  correspcmding  measurement-prediction  comparisons  are  presented  in 
appendix  B  (pages  B-22  to  B-33).  Very  few  of  the  lONCAP  and  none  of  the  ASAPS  predictions 
demonstrated  two-hop  propagation  at  this  range. 


Figure  5.  Twelve-Month  Frequency  Availability  for  Cape  Schmidt  and  Loparskaya 


CONCLUSIONS  AND  RECOMMENDATIONS 


This  rqxnt  has  summarized  woric  peifcvmed  to  extend  the  extrapolation  of  ionosonde 
measurements  from  the  single-hop  oblique  geometry  to  the  two-hop  case.  A  computer  program 
called  IONOLNK2  was  develq)ed  frcra  the  one-hop  program  lONOLINK  to  permit  FOT,  MUF, 
and  HPF  con^arisons  widi  lONCAP  and  ASAPS  predictions. 

The  IONOLNK2  program  should  be  used  for  mid-ladtude  two-hop  paths  in  which  more 
simultaneous  data  are  available  and  fewer  E-layer  reflections  are  observed.  This  effort  should  be 
performed  with  both  the  lONOLINK  and  IONOLNK2  programs  to  provide  a  complete 
characterization  of  frequency  availability  for  the  one-  and  two-hop  propagation  paths,  respectively. 
Once  the  investigation  is  ccxnpleted,  and  the  relationship  between  empirical  extrapolation  and 
computer  prediction  is  firmly  established  fra*  mid-latitude  paths,  an  extensive  high-latitude  study 
enconqMissing  more  years  of  ionosonde  data  should  be  performed.  This  high-latimde  effort  should 
include  use  of  the  ICECAP  computer  model,  so  that  a  comparison  of  its  predictions  with  one-  and 
two-hq)  extrapolatitm  as  well  as  with  lONCAP/ASAPS  predictitms  may  be  made. 
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APPENDIX  A 

I0N0LNK2  COMPUTER  PROGRAM 


A.1  DESCRIPnON  AND  INSTALLATION 

IONOLNK2  is  a  Microsoft  FORTRAN  S.O  computer  program  developed  and  executed  on  a 
Gateway  2000  486/33-MHz  PC.  The  program  was  written  and  dd)ugged  using  the  Power  Workbench 
^WB)  development  tool  from  Microsoft  (MS)  which  is  available  with  its  FORTRAN  product  MS- 
FORTRAN  5.0  employs  ANSI  FORTRAN  77  as  well  as  several  additional  features  for  compatibility 
with  VAX  FORTRAN  avulable  from  the  Digital  Equipment  Corporation.  The  purpose  of  the 
IONOLNK2  program  is  to  accumulate  numerical  densities  for  two-hop  oblique  standard  maximum 
usable  frequency  (MUF  )  values  using  input  files  in  the  "NEW  URSI"  data  format  as  published  in 
Intemationid  Council  of  Scioitific  Unions  Panel  on  World  Data  Centers  "GUIDE  to  &e  WORLD 
DATA  CENTER  SYSTEM",  Part  2.  These  parameters  include  hourly  junction  frequencies  (JFs), 
"standard"  maximum  usable  frequency  (MUF)  multipliers,  and  ionospheric  reflection  heights  (virtual 
heists),  for  the  each  of  the  observed  ionospheric  layers. 

IONOLNK2  consists  of  a  two  source  code  modules,  IONOLNK2.FOR  and  HFLINK2.FOR  which 
are  compiled  and  linked  by  the  PWB  tool  using  tbe  "make"  file  called  IONOLNK2.MAK  into 
IONOLNK2EXE.  Proper  use  of  the  MB-PWB  tool  is  beyond  the  scope  of  this  iqipendix.  The 
executable  image  is  invoked  at  the  user  prompt  from  drive  C  by  the  command: 

OIQHQLNK? 

where  the  highlighted  and  underlined  text  indicates  user  entry.  The  ”0"  in  this  example  is  displayed 
by  the  PC.  The  program  dien  asks  the  user  for  the  name  of  the  desired  IONOLNK2  input  file  with  the 
following  prompt: 

Type  name  of  input  file  widiout  extension  ".INF” 

The  user  then  types  frie  name  of  the  iiqiut  file  MiAout  die  trailing  (.)  or  file  extensioiL  The  IONOLNK2 
program  assumes  that  the  input  file  extension  is  ".INP".  For  example,  if  the  user's  input  file  is  called 
DKRS1991.INP  and  it  resides  in  directory  "INPUT"  on  drive  "D”,  then  the  user  would  respond  to  the 
prompt  with 

and  die  IONOLNK2  program  would  automatically  read  file  DKRS1991.INP  in  the  indicated  directory. 

AJl  INPUT  FILE  DESCRIPTION 
A^l  IONOLNK2  Input  Fde 

The  user-input  file  must  be  created  using  a  standard  ASCII  editor.  IONOLNK2  reads  this  input  file 
in  a  format-fi’ee  format,  so  exacdy  one  comma  (,)  or  one  or  more  spaces  may  be  used  to  delimit  entries. 
For  example,  consider  the  following  example  input  file: 


•DnCSON*,  •RESOLUTE  BAY . 1 

1. 1, 1, 1 . 2 

24,31,14,2000  . 3 

T,T,T,T,T,T,T,T,T,T,T,T,T,T,T,T,T,T,T,T,T,T,T,T,T  ....  4 

T,T,T,T,T,T,T,T,T,T,T,T . 5 

F . 6 

T,T,T,T,T,T,T,T,T,T,T,T,T,T,T,T,T,T.T,T.T.T,T,T,T,T,T,T,T,T,T  7 

T,T,T,T,T,T,T,T,T,T,T,T,T,T,T,T,T,T,T,T,T,T,T,T  ....  8 

10,20.0,40.0,60.0,80.0,100.0,200.0.400.0,800.0,1000.0,2000.0  9 

2.0, 32.0, 60  .  10 

100.0,4000.0,60  .  11 

0.0, 90.0,  60  .  12 

1, T.  0.0, 12.0,60 . 13 

2,  F,  0.0, 30.0, 60  .  14 

3,  F,  0.0, 30.0, 60  .  15 

4,  T,  2.0, 4.0, 50 . 16 

5,  F,  0.0, 600.0,  60  .  17 

6,  F,  0.0, 600.0, 60  .  18 

7,7,0.0,60.0,60  .  19 

8,  F,  0.0, 30.0, 60  .........  20 

9,  F,  0.0, 10.0,50 . 21 

10,  F,  0.0, 600.0,  edf . 22 

11,1,150.0,600.0,45  .  23 

12,1,0.0,5.0,25 . 24 

13,  F.  1.0, 30.0, 60  .  25 

14,  T,  90.0, 150.0,  60  .  26 

15,  F,  0.0, 600.0,  60  .  27 

16,  T,  0.0, 10.0,  50 . 28 

17,  F,  1.0, 30.0, 60  .  29 

18,7,0.0,5.0,25 . 30 

19,  F,  1.0,30.0,60  .  31 

20, 7,  90.0, 190.0, 20 .  32 

21,  F,  1.0, 30.0,  60  .  33 

22, 7, 0.0, 5.0, 25  .  34 

23,  F,  0.0, 10.0,  50 . 35 

24,  F,  0.0, 600.0, 60  .  36 

25,  F,  1.0,30.0,60  .  37 

26,7,0.0,16.0,40  .  38 

27,  F,  1.0,30.0,60  .  39 

28,  F,  1.0,30.0,60  .  40 

29,  F,  1.0,30.0,60  .  41 

30,  F,  1.0,30.0,60  .  42 

F . 43 

’D;VNUWaON\DK37391.NEW’, 'D:\NUWCION\RS9749I.NEW’  44 

record  (line)  is  shown  numbered  to  the  right  of  the  input  values  to  provide  an  index  for  reference 
to  eadh  item  in  diis  description.  These  numbers  must  not  i^pear  in  the  user-input  file  or  an  input  error 
vnll  occur.  A  description  of  each  input  value  referenced  by  these  numbers  is  provided  in  Table  A-1 . 
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Table  A.1  Input  Variable  Descriptions  for  10N0LNK2  Program 


INPUT 

FILE 

ROW 

NUMBER 

ITEM 
NUMBER 
IN  ROW 

FORTRAN 
VARIABLE  TVre 

DESCRIPTION 

TYPICAL 

VALUES 

1 

1.2 

CHARACTER*20 

Names  of  two  ionosonde  stations 
used  to  find  the  start  of  each 
month's  data  records 

mSSm 

2 

1 

INTEGERS 

First  hour  in  time  interval  for  bin- 
accumulation 

1 

m 

2 

M 

First  day  in  time  interval  for  bin 
accumulation 

1 

m 

3 

« 

First  month  in  time  interval  for 
bin  accumulation 

1 

" 

4 

« 

First  year  in  time  interval  for  bin- 
accumulation 

1970 

3 

1 

INTEGERM 

Last  hour  in  time  interval  for  bin- 
accumulation 

21 

m 

2 

11 

Last  day  in  time  interval  for  bin 
accumulation 

31 

m 

3 

n 

Last  month  in  time  interval  for  bin 
accumulation 

12 

m 

4 

N 

Last  year  in  time  interval  for  bin- 
accumulation 

4 

1-25 

LOGICAL*! 

Use^gnore  MASK  for  up  to  25 
years  of  ionosonde  data  files,  one 
file  per  year,  e.g.,  1970  -  1972 

lorl 
for  each 
entry 

5 

1-12 

LOGICAL*  1 

Use/ignore  MASK  for  up  to  12 
monAs  of  ionosonde  data  for  each 
year  (file),  1  tiirough  12 

lor£ 
for  each 
entry 

6 

1 

LOGICAL*! 

Logic  switch  for  program 
IONOLNK2,  not  used  in  program 
IONOLNK2;  T  creates  hourly 
FOT,  MUF  (standard),  and  HPF 
stafistics  for  oblique  HF  links,  F 
oeates  same  density  tabulations 

lorE 

ignored 

7 

1-31 

LOGICAL*! 

Use^gnore  MASK  for  up  to  31 
days  of  ionosonde  data  for  eadi 
month  (1  through  3 1) 

lorE 
for  each 
entry 

8 

1-24 

LOGICAL*! 

Usefignore  MASK  for  up  to  24 
hours  of  ionosonde  data  for  each 
day  (1  through  24) 

lorE 

for  each 
entry 

.-3 


Table  A.1  Input  Variable  Descriptions  for  10N0LNK2  Program  (Continued) 


INPUT 

FILE 

ROW 

NUMBER 

riEM 
NUMBER 
IN  ROW 

FORTRAN 
VARIABLE  TYPE 

DESCRIPTION 

VALUES 

9 

1 

INTEGER*4 

Input  to  track  lONOLINK  []  input 
files  but  not  used  in  IONOLNK2 
(10  MAX) 

10 

9 

2-11 

REALM 

Oblique  path  lengths  (km)  for 
lONOLINK  [],  mi  used  in 
program  IONOLNK2  (10  values 
MAX) 

&0 

to 

4000.0 

10 

1 

REALM 

Lower  frequency  limit  (MHz)  for 
IONOLNK2-computed  FOT, 

MUF,  &  HPF  values 

2^ 

M 

2 

REAL*4 

Upper  frequency  limit  (MHz)  for 
IONOLNK2-computed  FOT, 

MUF,  &  HPF  values 

60.0 

N 

3 

INTEGERM 

Number  of  frequency  bins  for 
IONOLNK2-computed  FOT, 

MUF,  ft  HPF  values  (60  ^iAX) 

H 

4 

LOGICAL*  1 

Lo^cal  switch  to  turn  off  sporadic 
£  layer  in  lONOLINK  Q,  mt  used 
inIONOLNK2 

lorE 

11 

1 

REALM 

Lower  height  limit  (km)  for 
IONOLNK2-compiited  FOT, 

MUF,  ft  HPF  values  in 
lONOLINK  Q,  mt  used  in 
IONOLNK2 

M 

N 

2 

REALM 

Upper  height  limit  (MHz)  for 
IONOLNK2-computed  FOT, 

MUF,  ft  HPF  values  in 
lONOLINK  Q,  mt  used  in 
IONOLNK2 

600.0 

N 

3 

INTEGERM 

Number  of  height  bins  for 
IONOLNK2-computed  FOT, 

MUF,  ft  HPF  values  (60  MAX)  in 
lONOLINK  Q,  mt  used  in 
IONOLNK2 

00 

Table  A.1  Input  Variable  Descriptions  for  IONOLNK2  Program  (Continued) 


INPUT 

FSX 

ROW 

NUMBER 

ITEM 
NUMBER 
IN  ROW 

FORTRAN 
VARIABLE  TYPE 

DESCRIPTION 

VALUES 

12 

1 

REALM 

Lower  antenna  elevation  angle 
limit  (deg)  for  IONOLNK2- 
comput^  FOT,  MUF,  &  HPF 
values 

M 

m 

2 

REALM 

Upper  antenna  elevation  angle 
limit  (deg)  for  IONOLNK2- 
comput^  FOT,  MUF,  &  HPF 
values 

90.0 

M 

3 

INTEGERM 

Number  of  angle  bins  for 
IONOLNK2-computed  FOT, 

MUF,  &  HPF  values  (60  MAX) 

60 

13. 14. 
15. 19. 
20. 24. 
25. 28, 
29, 30, 
31,33, 
34, 37, 
38, 39 

1 

LOGICAL*  1 

Logical  switch  to  turn  on  analysis 
for  the  corresponding  ionosonde- 
scaled-frequency  parameter  codes 
described  in  Section  A.2.2.  A  "T” 
tells  the  IONOLNK2  program  to 
process  the  parameter,  not 
processed  if  "F". 

Ior£ 

M 

2 

REAL*4 

Minimum  value  of  frequency  to  be 
used  in  density  accumulation  bins 

M 

n 

3 

REAL*4 

Maximum  value  of  frequency  to  be 
used  in  density  accumidation  bins 

30.0 

n 

4 

INTEGER*4 

Number  of  firequency  bins  to  be 
used  in  density  accumulation  bins 

u 

16. 21, 
35,40 

1 

LOGICAL*! 

Logical  switch  to  turn  on  analysis 
for  the  corresponding  ionosonde- 
scaled-factor  parameter  codes 
described  in  Section  A.2.2.  A  ”T" 
tells  the  IONOLNK2  program  to 
process  the  parameter,  not 
processed  if  "F". 

1 

lorE 

N 

2 

REAL*4 

Minimum  value  of  factor  to  be 
used  in  density  accumulation  bins 

M 

N 

3 

REAL*4 

Maximum  value  of  factor  to  be 
used  in  density  accumulation  bins 

Hi 

N 

4 

INTEGER*4 

Number  of  factor  bins  to  be  used  in 
density  accumulation  bins 

20 

Table  A.1  Input  Variable  Descriptions  for  10N0LNK2  Program  (Continued) 


INPUT 

FILE 

ROW 

NUMBER 

ITEM 
NUMBER 
IN  ROW 

FORTRAN 
VARIABLE  TYPE 

DESCRIPTION 

VALUES 

17, 18. 
22, 23, 
26, 27, 
32, 36 

1 

LOGICAL*  1 

Logical  switch  to  turn  on  analysis 
for  the  corresponding  ionosonde- 
scaled-height  parameter  codes 
described  in  Section  A.2.2.  A  "T" 
tells  the  IONOLNK2  program  u> 
process  the  parameter,  not 
processed  if  "F". 

lorZ 

M 

2 

REALM 

Minimum  value  of  height  to  be 
used  in  dmsity  accumulation  bins 

M 

N 

3 

REALM 

Maximum  value  of  height  to  be 
used  in  density  accumulation  bins 

600,0 

II 

4 

INTEGER*4 

Number  of  height  bins  to  be  used 
b  density  accumulation  bins 

60 

41,42 

1 

LOGICAL*! 

Lo^cal  switch  to  turn  on  analysis 
for  the  coirespondmg  ionosonde- 
scaled  parameter  code  as  described 
b  Section  A.2.2.  A  T"  tdls  the 
IONOLNK2  program  to  process 
the  parameter,  not  processed  if  *F*. 

lorE 

m 

2 

REAL*4 

Minimum  value  of  parameter  to  be 
used  b  density  accumulation  bins 

M(?) 

n 

3 

REAL*4 

Maximum  value  of  parameter  to  be 
used  b  density  accumulation  bins 

mo(?) 

m 

4 

INTEGER*4 

Number  of  parameter  bins  to  be 
used  b  density  accumulation  bins 

60(?) 

43 

1 

LOGICAL*! 

Logical  switch  to  have  densities  on 
total  elapsed  time  or.ek^sed  time 
for  which  ionosonde  data  was 
possible,  i.e.,  no  equipment  failure 
(qualifybg  letter  "C") 

Ior£ 

44-68 

1 

CHARACTER*20 

Ionosonde  data  filename  pairs, 
bcb^g  drive  and  subdirectory, 
wib  up  to  25  filename  pairs 
allowed 

See 

Section 
A.2.1  for 
examples 
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Table  A.2  lonosonde  (New  URSI)  Data  Format 


Record 

Columns 

Description 

1 

1-20 

Station  name 

1 

21-25 

Station  code 

1 

26-29 

Standard  time  meridian  of  the  station  (e.g. 
150W.  90E,  etc.,  with  OOOW  or  OOOE  =  UT) 

1 

30-  33 

Geographic  co-latitude  in  t«iths  of  a  degree 

1 

34-37 

Geognqihic  East  longitude  in  tenths  of  a  degree 

1 

38-41 

Year 

1 

42-43 

Month 

1 

44-45 

Parameter  code 

(See  UAG-23  characteristic  code) 

1 

46-120 

Spare 

2 

24  X  5-character  code 

Hourly  data  for  the  first  day  of  the  month  (i.e. 
foF2  values  -  078  R,  080  ,  etc.) 

3 

24  X  5-character  code 

Hourly  data  for  the  second  day  of  the  month 
(i.e.  foF2  values  -  078  R,  080  ,  etc.) 

24  X  5-character  code 

• 

24  X  5-character  code 

, 

24  X  5-character  code 

» 

32 

24  X  5-character  code 
(if  available) 

Hourly  data  for  the  thirty-first  day  of  the  month 
(If  less  than  3 1  days,  blank  fill.) 

33 

Medians 

34 

Median  Count 

35 

Upper  Quartile 

36 

Lower  Quartile 

37 

Upper  Decile 

38 

Range 

39 

Lower  Decile 

40 

Spare 

Not  all  of  the  five<haracter  scaled-parameter  values  woo  reported  for  eadi  hour  of  the  ionosonde 
da*"  provided  for  the  study.  In  this  case,  blank  fill  was  used  as  a  place  holder.  Note  that  the  files  used 
in  this  study  were  morged  from  separate  parameter-mondi  files  to  form  complete  parameter  files  for  all 
available  months.  The  completed  ionosonde  data  files  were  named  CS66991.NEW,  DK37391.NEW, 
LY168^1.NEW,  and  RS97491.NEW  for  Cape  Sdimidt,  Dikson,  and  Loparskaya  in  Russia  and  Resolute 
Bay  in  Canada,  respectively.  None  of  these  files  included  rows  33-40  in  the  monthly  data  blocks  for 
each  parameter.  If  N  parameters  have  been  scal^  then  tiiere  will  be  N  blocks  for  each  of  the  months 
mntainwH  within  a  single  file.  Also,  each  pair  of  files  used  in  the  two-hop  IONOLNK2  analysis  must 
ffnntain  exactly  the  same  parameter-month  blocks  or  program  execution  will  terminate.  Thus,  "null" 
blodcs  may  have  to  be  inserted  into  files  that  include  the  fust  record  (see  Table  A.2)  followed  by  31 
rows  of  24  "equipment  fault"  entries  given  by  "  C  ".  For  the  standard  N  -  14  parameters,  the 

ionosonde  file  should  be  organized  as  shown  in  Table  A.3.  A  complete  description  of  the 
codcs,  desciiptive,  and  qualification  letters  summarized  in  Table  A.4a,  b,  and  c,  respectively, 
can  be  found  in  UAG-23. 
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Table  A3  lonosonde  Data  File  Organization 


Block 

Month 

Parameter 

(Charaaeristic) 

Numeric 

code 

1 

January 

foF2 

00 

2 

January 

M(3000)F2 

03 

13 

January 

fxl 

51 

14 

January 

fml 

52 

15 

Fd)ruary 

foF2 

00 

16 

February 

M(3000)F2 

03 

167 

December 

fxl 

51 

168 

December 

Bnl 

52 

Table  A.4a  Parameter  Codes 


F2 


F 


E  layer 


Es 


Outer 


Spread 


TEC 


00-f6F2 

01-ficF2 

02-feF2 

10-foFl 

11-fidFl 

13-M(3000)F1 

20-foE 

22-foE2 

24-hE 

30-foEs 

31-fxEs 

32-fbEs 

40-foF1.5 

42-finin 

43-m(3000)F1.5 

50-foI 

5I-ficI 

52-finI 

70-1(2000) 

71-1 

72-I(xxx) 

14-hTl 


26-h'E2 


33-fEs 


44-hT1.5 


Table  A.4b  Quailing  Letters 

Letter _ Meaning _ 

A  Less  than  (used  only  in  case  of  total  blanketing) 

D  Greater  than _ 

E  Less  than _ 

I  Interpolated _ 

J  Deduced  from  x  component _ _ 

M  Mode  uncertain _ 

O  Deduced  from  o  component _ 

T  Smoothed  bom  sequence _ 

U  Uncertain _ 

Z  Deduced  from  z  component _ 


Table  A.4c  Descriptive  Letters 


mnrm 

Meaning 

A 

Blanketing 

B 

Absorption 

C 

Non-ionospheric  (equipment) 

D 

Above  upper  freq.  range 

E 

Below  Iowa  freq.  range 

F 

Spread  echoes 

G 

Ionization  density  too  small 

H 

Stratification 

K 

Night  E  laya  present 

L 

Insufficiently  defmed  cusp 

M 

Mode  uncertain 

N 

Superimposed  layers 

0 

Measurement  refers  to  o  component 

Q 

Range  spread 

R 

Attenuation  nea  critical  freq. 

S 

Interfaence 

T 

Interpolated 

V 

Forked  trace 

W 

Above  height  range 

X 

Measurement  refers  to  x  component 

Y 

Lacuna  (tilt) 

Z 

Measurement  refers  to  z  component 

AJ3  OUTPUT  FILE  DESCRIPnON 

The  IONOLNK2  computer  prognun  can  genoate  two  different  output  file  formats  depending  on  the 
value  of  input  file  record  number  six  (6).  In  either  case,  the  ou^t  file  begins  with  a  description  of  the 
two-hop  link  geome&y,  including  the  geographic  coordinates  of  the  two  ionosonde  stations.  If  record 
six  is  "false"  ("F"),  the  IONOLNK2  output  filename  has  the  same  filename  as  the  input  file  but  with  the 
extension  ".OUT"  instead  of  ".IMP".  It  contuns  a  tabular  normalized  number  density  (probability 
density  function)  for  each  possible  propagation  mode  (currently  one  of  the  16  possible  two-hop  modes 
F2-F2,  F2-F1,  F2-E,  F2-Es,  F1-F2,  Fl-Fl,  Fl-E,  Fl-Es,  E-F2,  E-Fl,  E-E,  E-Es,  Es-F2,  Es-Fl,  Es-E,  or 
Es-Es),  the  average  standard  MUF  v«iue,  ;dnd  tjie  standard  deviation  of  the  standard  MUF  value  for  each 
hmir  of  the  day.  These  statistics  are  created  ova-  the  entire  time  interval  requested  by  the  input  file. 
The  format  for  these  density  tables  is  depicted  in  Tables  A.Sa  and  A.5b.  A  sample  output  file  called 
DKRS1991.0UT  (Dikson-Resohite  Bay)  is  provided  with  the  software  package. 

The  second  ouqiut  file  format  was  developed  to  provide  a  direct  output  for  comparison  with 
lONCAP,  ASAPS,  or  otha  HF  prediction  code  ouq>ut  This  output  file  provides  the  optimum  working 
fir^ency  (OWF  or  FOT,  10%  not  exceeded),  the  median  standad  MUF  value  (S0%  not  exceeded),  and 
tile  Iowa  decile  standad  MUF  value  (90%  not  exceeded)  versus  hour  of  day.  In  addition,  the 
corresponding  average  antenna  take-off  angle  for  each  standad  MUF-distribution  value  is  also 
provided.  Thi.q  output  file  has  the  extension  ".LNK"  and  cannot  be  aeated  simultaneously  with  the 
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corresponding  ".OUT*  file.  Parameter  six  (6)  in  the  input  file  determines  the  form  of  IONOLNK2 
ou^Hit.  The  format  of  this  file  is  outlined  in  Table  6  and  a  sample  "".LNK"  file  is  provided  with  the 
software. 


Table  A.5a  Propagation  Mode  Density  Table 


Hour  Intervals  (UT) 

PROPAGATION 
MODE:  ... 

0-1 

1-2 

2-3 

23-24 

NET 

100  Pr 

wrmm 

nm 

nnn 

nnn 

nnn 

nnn 

WffSm 

AVGFRQ 

fff 

ff 

fff 

fff 

fff 

fff 

SDVFRQ 

■ii 

sss 

sss 

sss 

sss 

sss 

sss 

sss 

^nwi  is  a  probability  expressed  as  a  percentage  (0  ^  mm  ^  ICO) 


a  frequency  (MHz)  ^sss  is  frequency  standard  deviation  (MHz) 


Table  A.5b  Standard-MUF  Density  Table 


Hour  Intervals  (UT) 

Frequency  (MHz) 

0-1 

1-2 

2-3 

• 

. 

, 

23-24 

NET 

ram 

ram 

ram 

ram 

ram 

ram 

ram 

ram 

fiA 

ram 

root 

ram 

ram 

ram 

ram 

ram 

ram 

fyfA 

ram 

ram 

ram 

ram 

ram 

ram 

mm 

ram 

. 

ram 

ram 

ram 

ram 

ram 

ram 

ram 

ram 

. 

ram 

ram 

ram 

ram 

ram 

ram 

ram 

ram 

. 

ram 

ram 

ram 

ram 

ram 

ram 

ram 

ram 

ram 

ram 

ram 

ram 

ram 

ram 

ram 

ram 

indicate  probability. 


Table  A.6  FOT,  MUF,  HPF  Table  in  "^XNK"  Output  File  for  Each  Month 


10%  NOT  EXCEEDED  (FOT) 

year 

month 

(HOUR;) 

0-1 

1-2 

2-3 

• 

• 

• 

23-24 

NET 

FOT: 

ff 

fff 

M 

fff 

ff 

ff 

M 

ANG: 

OOtP’ 

OQQ 

Qcta 

aaci 

aaa 

aaa 

aaa 

aaa 

50%  NOT  EXCEEDED  (MUF) 

year 

month 

1 

(HOUR;) 

0-1 

1-2 

2-3 

• 

. 

. 

23-24 

NET 

MUF; 

fff 

fff 

ff 

fff 

ff 

m 

fff 

ANG: 

QCUZ 

acta 

QCUJ 

aaa 

aaa 

aaa 

aaa 

aaa 

1  90%  NOT  EXCEEDED  (HPF) 

year 

month 

(HOUR;) 

0-1 

1-2 

2-3 

. 

- 

. 

23-24 

NET 

HPF. 

ff 

Jff 

iff 

ff 

Jff 

fff 

fff 

fff 

ANG; 

cuux 

aaa 

aaa 

aaa 

aaa 

aaa 

aaa 

W*  frequency  in  MHz  '^aaa  -  angle  in  degrees 
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A3  IONOLNK2  PROGRAM  USTING 

A3.1  10N0LNK2.F0R  Source  Program 

SFREEFORM 

SLARGE 

INTERFACE  TO  SUBROUTINE  HFLINK2(  il,  i2, 11. 12  ) 
INTEGERS  il,  i2 
LOGICAL*!  11(4,4),  12 
END 

INTERFACE  TO  SUBROUTINE  PARMBN2(  il,  i2.  i3,  i4,  rl,  cl,  c2  ) 
REAL*4  rl 

INTEGER*4  il.  i2.  i3,  i4 
CHARACTER*!  cl,  c2 
END 


"  Program  HFSOND2 

N 

*  Purpose:  Process  an  lONOSONDE  data  file  following  the  conventions 

”  established  in  UAG-23  and  UAG-23A  and  detennine  the  optimum 
”  frequencies  for  an  HF  communications  link  with  its  midpa^ 

"  above  the  ionosonde. 

•I 

”  Assumptions:  It  is  assumed  that  both  ionosonde  data  files  have  the 
"  same  number  of  parameter  data  blocks  per  month.  If  not, 

”  then  these  input  files  must  be  modified  to  include  a 

*  "filler"  block  containing "  C "  entries  for  the 

"  3 1  rows  of  24  scaled-parameter  values. 

PROGRAM  HFSOND2 

"  Declare  variable  types  and  dimension  arrays 

"  Define  s.  general  input  record  with  120  columns  and  a  bolder  for  the 
"  previous  record 

CHARACTER*120  RECORD J(2),  RECORD(2) 

"  Define  parametms  to  hold  values  describing  the  ionosonde  station 

CHARACTER*20  STATIONjiame(2),  STATION_check(2) 
CHARACTER*5  STATION_code(2) 

CHARACTER*4  TIME  meridian(2),  STATION_COLAT_ddg(2) 
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CHARACTERS  STATION  ELNG  ddg(2).  YEAR(2) 

CHARACTER*2  MONTH(2),  OLD  MONTH 
CHARACTERS  f  RANGE(12) 

INTEGERS  P  codc(2),  N  file,  M  1,  M  2.  Y  1,  Y  2.  Y_3.  Y  4 
INTEGER*4  YEAR  i,  MOimi  i,  DAY  i.  HOUR  i 
INTEGER*4  YEAR J,  MONTH J.  DAY J.  HOURJ 

*  Define  S-character  code  derived  from  UAG-23  for  hourly  interval 

CHARACTER*  1  P_char(3),  Q_char,  D_char 

"  Define  bins  for  yearly  distribution  of  hourly  values 

REAL*4d_P,P  l.P  2,P  3,P  value,  MAX  hours(l;30,l;24) 

REAL*4  PJowO:30X  P_hgh(l730).  P_scale(i:3) 

INTEGER*4i_FILE.j  CODE(2),N  CODE,  i  CODE,P  index(l;73),N_DISP 
INTEGER*4N  Pbins(T:30),i  MONmPARM_code(2)7P_type(l:30) 
INTEGER*4  ERR_LINE(2) 

"  Declare  arrays  for  HF  link  calculations 

REAL*4  fMHzJow,  fMHzJigh,  dkmjow,  dkmjhgh,  Rejbn,  Rearfrijon 

REAL*4  fMHzJhr^^p,  fMHz_disp 

REAL*4  ad^low,  adg^hgh,  d JMHz,  djlkm 

REAL*4  ajimit,  adgjbin,  adgjhrjbin,  adgjnodejbin 

REAL*4  adgjnodejirjbin,  AVGjuIgjtnoideJir,  AVGjadgjnode 

REAL*4  d_iiinit,  dkmjbin,  dkm^hrjbin,  dkm_modeJbin,  D_km(10) 

REAL*4  dkm_mode_hrJbin,  AVG_dkm_mode_hr,  AVG_dkmjnode 

REAL*4  f_limit,  fMHzJbin,  fMHzJurJbin,  fMHz_modeJbin 

REAL*4  fMHzjmode_hrJbin,  AVG_fMHz_mode_hr,  AVG_fMHz_mode 

REAL*4  SDV_fMHz_inode_hr,  SDy_fMHz_mode,  N_noLOS 

REAL*4  N_dy,  NJir_dy,  N_m^e_dy,  N_modc_hr_dy,  N_dy_hr,  N_hr,  N_GOOD 

REAL*4  COLAT_dg(2),  LNG_dg(2),  Ds_kin,  Asjdg,  azimuthjdg(2) 

INTEGER*4  i_DAY,  i_HOUR,  i  DST,  i  FRQ,  i  MODE,  i  SITE 
INTEGER*4  N  FRQ,  N  ANG,  N_DST,  N_Dkm 
INTEGER*4  N^bins,  month_OLD 

CHARACTER*5  L_heading(3),  MODE_name(l:4,l:4) 

CHARACTER*25  DIST_headmg 

"  Define  storage  for  each  parameter  code,  descriptive,  and  qualifying 
"  letter. 

REAL*4  foF2,  fxF2,  f2F2,  M3000F2,  hF2,  hpf2,  foFl,  fxFl 
REAL*4  M3000F1,  hFl,  hF,  foE,  foE2,  hE,  hE2,  foEs 
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REALM  fxEs,  fbEs,  fEs.  hEs.  foFldS,  finin,  M3000FldS,  hfldS 
REAL*4  fol,  fid,  fml,  12000,  H,  Ixxx 

Define  unused  logical  parameter  employed  in  the  lONSTATS  program 

LOGICAL*  1  DIV_flag 

Define  a  character  string  of  BLANK  Jine  for  null  string  searches 

CHARACTER*21  BLANK  Jine 

Define  list  of  input  files  for  analysis 

CHARACTER*25  SI  files(l:25),  S2  files(l:25),  INP_file 
CHARACTER*25  LNK.file,  FILE_NAME,  OUT_file,  file_INP 

Define  heading  variables  for  ou^ut 

CHARACTER*  10  P_names(l  :30) 

CHARACTER*S  QUAL  heading,  DESCJieading 
CHARACTER*5P  heading(l:3),DIM  headtng(l;4) 

CHARACTER*25  HOUR  heading 
CHARACTER*!  QUAL_value(l:10),  DESC^value(l:23) 

LOGICAL*!  YEAR_mask(!:25),MONTHjnask(!:!2),DAY  mask(l:3!) 
LOGICAL*!  HOUR  mask(!:24),  CODE  mask(l:30),  C0DE3ag 
LOGICAL*!  P  flag,  LINK  flag(4,4),  L  flag 
LOGICAL*!  pToJFLAG,  M0_FLAG,  ^_FLAG 
LOGICAL*!  RE^Jlag,  RUNJlag,  Reptjlag,  DAYJlag,  HOUR_flag 
LOGICAL*!  MONTHLY Jlag,  FIRSTjnonth,  GEOGJlag 

**«*4i****4i«****************«**********************4i*******4i**4i**4i4i4> 

Define  COMMON  blocks 

COMMON/IONSNP/foF2(2,3!,24),  M3000F2(2,3!,24),  hF2(2,31,24),  - 
foF!(2,3!,24),  M3000F! (2,31,24),  hF!(2,31,24),  - 
hF(2,3 1,24),  foE(2,3 1,24),- 
hE(2,3!,24),  foEs(2,3!,24),  fbEs(2,3!,24),  - 
£Es(2,3!,24),  hEs(2,3!,24) 

fmin(2,3!,24),  fid(2,3!,24) 

COMMON  /MISC/  Re_km,  Rearth_km 

COMMON  /LIMIT/  fMHz Jow,  fMHz_hgh,  NJFRQ,  dkm Jow,  dkm_hgh,  N J>ST,  - 
adg  low,  adgjigh,  N_ANG,  LNK Jile,  dJMHz,  d_dkm,  - 
d julg,  As_rd,  Dsjcm 


*  COMMON  /ANGLE/  adg_bin(0:61),  - 

*  »dg^hr_bin(l;24,0:61),  adg_mode_bin(l:4,l;4,0:61),  - 

"  tdgjn^e  hr  bin(l:4,l:4,l:24,0;^l), - 

"  AVG  adgjnode  hr(l;4, 1:4, 1:24),  AVG_adg_mode(l:4,l  :4),  - 

*  ANGllirJrq(l:24,0:61,l:2) 

COMMON  /ANGLE/  ANG_hr_frq(l  :24,0:61,1 :2) 

"  COMMON  /RANGE/  dkm  bin(0:6I).  > 

"  dian_lir_bin(l:24,C:61).  dkin_^modc_bin(l:4, 1:4,0:61),  - 

"  dkm  m^e  hr_bin(l  :4, 1:4, 1:24, 0:61), - 

*  AVG_dkmjn^e_lir(l:4, 1:4, 1:24),  AVG_dkm_mode(l:4,l:4) 

COMMON /FRQBN/ f  limit(l:12,l:2),  fMHzJ)in(0:61),  - 

fMHz_hrJ)in(l:24,0:61),  fMHz_modeJbm(l:4, 1:4,0:61),  - 
fMHz  mode  hrjbiii(l:4,l:4,l:24,0:61),- 
AVG iMHz'mode  lir(l :4, 1 :4, 1 :24), - 
SDV  fMHz  mode  hr<l:4,l:4,l:24), - 

AVG  fMHz  mode(l:4,l:4),SDV_fMHz_mode(l:4,l:4),- 
fMHzl.lir_disp(l:24,l:12),  fMHz_disp{l:r2) 

COMMON /COUNT/ N  dy,N  hr  dy(l:24),N  mode  dy(l:4,l:4),N  GOOD,- 
N  mode  hr  dy(l:4,l-4,l:24),N  dy_hi(l:24),N  hr,- 
N^noLOS 

H  «*****4i*4i4i*********4i******«*«***********4i********4t***************** 

”  Initialize  constants 

DATA  P  index/l A3,4,5,6,4*0,7,8,0,9,10,0,1  1,3*0,- 

12,0,13,0,14,0,15,3*0,16,17,18,19,20,5*0,- 

21,0,22,23,24,5*0,25,26,27,17*0,28,29,30/ 


DATA P  type/3*l,3,2*2,2*l,3,2Al,l,2,2,4*l,2,- 
1,T,3,2,3*1,3*1/ 

DATA  P^scale/lO.0, 1.0, 100.0/ 

DATABLANK.line/'  V 

DATA  P  names/•foF2^•fxF2^T2F2^M(3000)F2^1lT2^•fp^^ 
•foFl','fitFr,’M(3000)Fl','h"Fr,'h"F,ToF,'foE2Vh"EVh"E2’.- 

tbEs^•fxEs^lbEs^m',•h"Es•,■foF1.5',•finin^'M(3000)F1.5•,- 

•h*fl.5','for,fxr,'finr,1(2000)’,T,T(xxx)V 

DATA  MODE  name/T2-F2',Tl-F2’,'  E-F2',*Es-F2’,  - 
f2-Fr,Tl-FlV  E-Fr,'Es-Fl',  - 
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•F2-EVF1-EV  E-EVEs-E',- 
T2-Es’,Tl-Es','  E-Es'.’Es-EsV 


DATA  P  heading/*  FRQ HGT  V 
DATA  DiM_hcadmg/'(MH2)’.'  (kin)*,TACTR‘,'(deg)V 
DATA  HOUR  heading  Unive^  Time  in  Hours  V 
DATA  QUALlieading/'  QUALV 
DATA  DESC.heading/*  DESC/ 

DATA  DIST  heading/Tath  Length  in  Kilometers'/ 

DATA  L^heatog/*  FRQ  •;  ANG HGT  Y 
DATAfJimit/2.0, 6.0, 10.0, 14.0, 18.0,  - 

22.0, 26.0, 30.0, 40.0, 60.0, 100.0, 150.0,  - 
6.0, 10.0,  14.0, 18.0, 22.0, 26.0, 30.0,  - 
40.0, 60.0,  100.0, 150.0, 300.0/ 

DATANJ)ISP/12/ 

DATA  f_RANGE  /*  2-6  '  6-10 '10-14 ',  '14-18 ',  *18-22 ',  - 

•22-26 ',  *26-30 ',  *30-40 ',  *40-60 ',  *60-100',  - 
to  150',  to  3007 

Reaith_km  *  6370.0 

Re  Jem  *  4.0  *  Reaithjm  /  3.0 

«  **«**4i*«4>*********«**«*«**««***«*«****4i********4i***********4>******» 


"  Open  ii^mt  file 

WR1TE(6,*)  Type  name  of  input  file  'whh  extension  ".INP"' 
READ(6,600)fileJNP 

"  Concatenate  input  and  ouQnit  file  names 

INPJUe  -  file_INP(l  :LEN  TRIM(file_INP))//'.INP 
OUT_file  =  file  INP(liEN  TRIM(file  INP))/AOUT 
LNK.file  «  fileIlNP(l  :LEN_TRIM(file_INP))//'.LNK' 

"  Read  input  file 

OPEN(Uim'=2,STATUS='OLD'JFILE*INPJile) 
READ(2,*)  STATION  check(l),  STATION_check(2) 
READ(2,*)HOUR  i,DAY  i,  MONTH^i,  YEARJ 
READ(2,*)  HOUR  f,  DAY_f,  MONTHJ,  YEARJ 
READ(2,*)  (  YEAR  maskO),  j  «  1, 25  ) 

READ(2,*)  (  MONTH  maskO),  j  »  1, 12  ) 

READ(2,*)  MONTHLY  flag 
READ(2,*) (DAY  maskO), j “  1, 31 ) 


READ(2,*)  (  HOUR  mask(j).  j  -  1, 24 ) 

READ(2,*)N  Dkm.  ( D  kmO),  j  *  1.  N  Dkm  ) 

R£AD(2.*)fI^  loy/.fMHz  hgh,N  FRQ 
READ(2,*)  dkm Jow,  dkm^h^  N.DST 
R£AD(2,*)  adgjow,  adg_^,  N  ANG 
DOi-1.30 

READ(2.*)  j.  CODE  maskO),  P  low(j).  P  hgh(j).  N_PbinsG) 
END  DO 

READ(2.*)DIV  flag 
"  READ(2.*)N_nLE 

"  Read  names  of  standard  ionosonde  files.  If  an  error  occurs  or  if 
"  if  no  input  filenames  are  present,  output  the  appropriate  message 
*  and  stop  execution 

LFBLE-l 

100  READ(2,*.ERR-998.EN1>999)  Sl_files(i_FILE).  S2_files(i_FILE) 

"  Determine  the  number  of  ionosonde  input  files 

DO  WHILE  ( i  FILE. GT.O) 
i  FILE«i  FILE+1 

MAD(2,*jENI>=102)  SI  files(i  FILE),  S2  files(i  FILE) 

END  DO 

102  N_file-i_FILE-l 
"  i_FILE-2 

"  DO  WHILE  (i  FILE  XE.N  FILE) 

"  i  FILE  =  i_Fn£+l 

"  READ(2,601)  Sl_files(i_FILE),  S2_fUes(i_FILE) 

"  ENDDO 

dJMHz  =  (  fMHz^hgh  -  fMHzJow  )  /  FLOAT(  N_FRQ  ) 
d  adg  =  ( adg_hgh  -  adg_low  )  /  FLOAT(  N_ANG  ) 
didkm  -  (  dkm Jigh  -  dkm  Jow  )  /  FLOAT(  N J>ST  ) 

"  INITIALIZE  all  bin  values  to  zero 

N_GOOD  =  0.0 
monthjold  *  0 
adgjbin  *  0.0 
adgjirjjin  *  0.0 
adg^modejiin  =  0.0 
adgjnodejirjiin  ®  0.0 
AVGjulgjmodeJhr  *  0.0 
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AVG_«dgjnodc  ■  0.0 
ANG_hr_frq  *  0.0 
dkmjbin  »  0.0 
dkmjhrjbin  “  0.0 
dkmjnodejbin  «  0.0 
dkmjmode_hrJbiii  •  0.0 
AVGjikm_mo<ic_hr  *  0.0 
AVG_dkin_inode  *  0.0 
fMHTbin-O.O 
fMHzJbrJbin  ■  0.0 
fMHzjmodeJbin  0.0 
fMHzjmodeJirJ>in  *  0.0 
AVG_fMHz_inodc_hr  «  0.0 
SDV_fMH2_modeJ»r  *  0.0 
AVG_£MH2_inode  “  0.0 
SDV_£MH2_modc  *  0.0 
fMHz^hrjdisp  *  0.0 
fMHz_disp  *  0.0 
N  dy*0.0 
NJir-O.O 
N_hrjdy  *  0.0 
N_dyjhr*0.0 
N_modc_dy  *  0.0 
N_modeJhr_dy  *  0.0 
N_months  «  0 

nIcode»o 


"  Initialize  all  ionosonde  parameter  storage  locations  to  >1.0.  Note  that 
”  the  commented  parameters  are  not  yet  used  in  the  calculation  of  forward 
"  propagating  signals. 


foF2»-1.0 
"  ficF2«-1.0 

"  feF2  =  -1.0 
M3000F2  *-1.0 
hF2  *-1.0 
"  hpf2*-1.0 

foFl*-1.0 
fidFl*-1.0 
M3000F1*-1.0 
hFl*-1.0 
hF*-1.0 
foE*-1.0 
"  foE2--1.0 

hE*-1.0 
"  hE2*-1.0 

foEs*-1.0 
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•  ibEs--1.0 
fE$--1.0 
I1ES--I.O 

•  fbFl(i5--1.0 

•  M3000Fld5--1.0 
"  hfldS-.l.O 

"  fol--1.0 

•  fid- -1.0 

"  finl--1.0 

"  12000- -1.0 

"  IXXX--1.0 


"  Print  heading  for  OUTPUT  file. 

m 

IF  (MONTHLY  flag)THEN 

OPEN(UNIT-8,STATUS-TJNKNOWN’^CCESS='APPEND’,FILE«LNK_file) 
WRITE(8,701)  LNK.file 
ELSE 

OPEN(UNIT-3,STATUS»'NEW,FILE-OUT  file) 

END  IF 

•  START  of  FILE  loop 

i  FUJE-O 
RUN^fiflg-.TRUE. 

C£OG  flag-EALSE. 

DO  WHILE  ( i_FILE  IT.  NJTLE  .AND.  RUN Jflag  ) 

”  Initialize  the  FIRST_month  flag  to  "TRUE",  indicating  that  the 
"  fiist  month  of  data  in  the  current  input  files  has  not  been  (reread 

"  AND  inaementtheinputfile  counter  to  a  "1". 

FIRST  month -.TRUE. 
i_FILE»'l_FILE+ 1 

"  Set  both  READ  Jflag  values  to  "T"  at  the  start  of  the  parameter  input  loop 
READJlag-.TRUE. 

"  Open  input  file  for  ionosonde  data  fiom  current  station 

OPEN(UNrr-4,STATUS='OLD',FILE-S  l_fiies(i_FILE)) 
OPEN(UNrr»7,STATUS='OLD',FILE=S2_files(i_FILE)) 


A-18 


Start  of  file-read  loop.  "R£AD_flag"  tells  the  program  to  keep 
reading  parameter  v^ues  for  the  current  YEAR,  MONTH,  and  DAY  in  the 
MONTH  if  set  to  T”.  If  "F",  then  it  discontinues  reading  the  current 
station  file. 


i  READ-0 
ERR  LINE(1)»0 
ERR  LlNE(2)-0 


DO  WHILE  (  READ_flag  .AND.  RUN_flag  ) 


i_READ  =  LREAD+l 

Read  the  next  line  of  the  current  input  file  and  check  the  first 
21  characters  to  make  sure  they  are  nonblank.  If  they  are  blank, 
then  a  new  record  should  be  r^. 

REC0RD_I(1X1:21)- BLANK  line 
DO  WHILE(  REC0RD_I(1X1:21)  BQ.  BLANK_lme  ) 
ERRJLINE(1)  -  ERR_LINE(1)  +  1 
READ(4,702,END-200)  RECORDJ(l) 

END  DO 

RECORD  1(2X1:21) -BLANK  line 
DO  WHIli(  RECORD J(2X1:21)  .EQ.  BLANK  line ) 
ERR_LINE(2)  -  E^_LINE(2)  +  1 
READ(7,702,END-200)  RECORDJ(2) 

END  DO 

Extract  station  name,  code,  standard  time  meridian,  geographic  co-latitude, 
east  longitude,  year,  month,  and  the  parameter  code 

DOi_SrrE»l,2 

STATION_name(i_SITE)  -  RECORDJ(i_SITEXl  :20) 


IF  (  STATION  name(i_SITE)  .NE.  - 
STATION  check(i  SITE)) THEN 
WRITE(6,*)' 

WR1TE(6,*) 'ERR:  Station  name  is  not  input  name!' 
WRITE(6,*)'l','  STATION  check:  ',  STATION_check(l) 
WRrrE(6,*) 'STATION  name:  ',  STATION_name(l) 
WRITE(6,*)  TILE  1  LINE: ',  ERR_LINE(1) 

WRITE(6,*) 

WRITE(6,*)  '2', '  STATION_check:  ',  STATION_check(2) 
WRITE(6,*)  'STATION  name:  ',  STATION_name(2) 
WRITE(6,*)  'FILE  2  LINE:  ’,  ERR_LINE(2) 


WRniE(6,*) '  Execution  has  been  halted.' 

WRrrE(6,*) '  Please  check  station  names  and  try  again.' 

STOP 

END  IF 

STATION  codc(i  SITE) -RECORD_I(i_SrreX2 1:26) 

TIME  meridian(i  SITE)  =  RECORDJ(i  SrrEX26:29) 
STATON  COLAT  ddg(i  SITE)  -  RECbRDJ(i_SITEX30;33) 
STATONlELNG_ddg(i_SITE)  »  RECORD J(i_SITEX34:37) 

END  DO 

*  Extract  the  YEAR  from  the  current  input  file.  If  the  YEAR  is  not  to 
"  be  included  in  the  analysis  (YEARjnask)  or  it  falls  outside  of  the 

*  desired  processing  interval,  then  exit  loop  which  reads  Ele  for  current 

*  YEAR  for  both  sites,  i.e.,  do  not  read  the  corresponding  file  for  the 

*  second  site. 

DOLSITE-1,2 

YEAR(i  SITE) -RECORD  I(i  SITEX38:41) 

Y  l-iaiAR(  RECORD  I(i  SrrEX38:38)  )  -  48 
y“’2  -  ICHAR(  RECORDlcf SrTEX39:39)  )  -  48 

Y  3  -  ICHAR(  RECORD  I(i  SniEX40:40)  )  -  48 
Yj4«ICHAR(  RECORD  10  SITEX41:41))-48 

i  YEAR- 1000 *Y  l  +  lob*Y2+10*Y3  +  Y4 
W  ( i  YEAR  XT.  YEAR  i  .OR  f  YEAR  .GT.  YE/AJ  OR- 
.NOT.  YEAR  maskO  FD-E)  )  THEN 
READ_flag-. FALSE. 

ENDIF 

END  DO 

IF  ( .NOT.  READJlag  )  CYCLE 

IF  (  YEAR(1)  .NE.  YEAR(2)  )  THEN 

WRITE(6,*)  "ERR  The  YEAR  values  differ  between  the  two’ 

WRnE(6,*) '  ionosonde  input  files!’ 

WRITE(6,'>) ’ACT:  Execution  has  been  halted.' 

WRnE(6,*)  Tlease  check  diese  files  and  try  again.' 

STOP 

ENDIF 

"  Extract  the  MONTH  from  the  current  input  file.  If  the  MONTH  is  not  to 
"  be  included  inthe  analysis  (MONTH_mask)  or  it  falls  outside  of  the 
"  desired  {xocessing  interval,  then  exit  loop  which  reads  file  for  current 
"  MONTH  for  both  sites,  i.e.,  do  not  read  the  corresponding  file  for  the 
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second  site.  If  this  is  the  fust  input  file,  initialize  MONTH  1  to 
i_MONTH. 

DOi  SITE- 1,2 

MONTHO  SITE) -RECORD  I(i  SrrEX42:43) 

M J  -  ICHAR(  RECORD  I(i  SrreX42:42)  )  -  48 
M,2  -  ICHAR(  record" I(rSITEX43:43)  )  -  48 
i  MONTH- 10  •M  1+M  2" 

^DO 

IF  (  MONTH(l)  .NE.  MONTH(2)  )  THEN 

WWTE(6,*)  'El^:  The  MONTH  values  differ  between  the  two' 

WRITE(6,*} '  ionosonde  input  files!' 

WRITE(6,*) 'ACT:  Execution  has  been  halted.' 

WRITE(6,*)  'Please  check  these  files  and  try  again.' 

STOP 
END  IF 

If  the  current  YEAR  is  the  fuud  year,  check  to  see  if  the  final  MONTH 
has  been  exceeded.  If  so,  stop  reading  parameter  blocks  for  the  current 
SITE  and  either  look  at  the  next  SITE'S  file  or  exit  to  process  the 
results. 

IF  ( i.YEAR  EQ.  YEAR  f  .AND.  i  MONTH  .GT.  MONTH  f  )  THEN 
READ  flag- JALSE. 

CYCLE 
END  IF 

If  the  current  MONTH  data  block  is  not  to  be  processed  or  if  die  first 
MONTH  data  block  to  be  processed  has  not  yet  been  encountered,  then  read 
the  next  31  records  to  skip  the  current  data  blodc  and  read  the  next 
data  block. 

ff  ( .NOT.  MONTH_mask(  i.MONTH  )  .OR.  - 

( i  YEAR  .EQ.  YEAR  i  .AND.  i  MONTH  .LT.  MONTHJ  )  )  THEN 
DO"i-l,31 

ERRJLINE(1)  -  ERR_LINE(1)  +  1 
READ(4,702,END=26b)  RECORD_I(l) 

ERR  LINE(2)  =  ERR  J.INE(2)  +  1 
READ(7,702END-26b)  RECORD_I(2) 

END  DO 
CYCLE 
END  IF 

IF  (  FIRST  month  )  THEN 
MONTHJ  -  i  MONTH 
FIRST  month- JALSE. 
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END  IF 


*  If  tills  u  the  first  file  read,  determine  the  latitude  and  longitude 

"  of  the  current  ionosonde  station  site  and  then  compute  the  inter-site 

*  Great  Circle  path  length  and  bearing  angles. 

IF  (  NOT.  GEOG.flag  )  THEN 
DOLSITE-1.2 

L  1-ICHAR(  RECORD  I(i_SITEX30;30)  )  -  48 
L  2-ICHAR(RECORD_I(i  SITEK31:31)).48 
l’3  -  ICHAR(  RECORD_I(rsmEK32:32)  )  -  48 
L_4-ICHAR( RECORD  I(i  SITEX33:33)).48 
I  DUM-100*L  1  +  10*L_2  +  L  3 
COLAT_dg(i_SnE)  -  FLOAT(  I_DUM  )  +  0.1  *  FLOAT(  L_4  ) 

L_1»ICHAR(  RECORD  I(i  SITEX34:34)  )  -  48 
L  2  -  ICHAR(  record’ I(i_SrrE)(35;35)  )  -  48 
L  3  *  ICHAR(  RECORDJ(i_SITEX36;36)  )  -  48 
L  4 -ICHAR(  RECORD  I(i  SITEX37:37)  )  -  48 
I  DIIM*=100*L  1  +  10’*L  2  +  L  3 
U^G_<*g(LSITE)  -  FLOAT(  F  DUM  )  +  0.1  •  FLOAT(  L_4  ) 

END  DO 

CALL  BEARING(  STAHONjiame,  COLAT_dg,  LNG_dg,  Dsjm,  - 
Asjrd,  azimuth_dg,  MONTHL YJlag  ) 


GEOG_flag  » .TRUE. 

END  IF 

"  Read  and  decode  the  parameter  code  whose  values  are  contained  in  the 
”  nesct-to-be-read  input  block  for  the  current  file  and  ionosonde  station. 

DO  i_SITE  *  1, 2 

P_code(l)  =  ICHAR(  RECORD_I(i_SITEX44:44)  )  -  48 
P_code(2)  =  ICHAR(  RECORD_I(i_SITEX45;45)  )  -  48 
PARM_code(i_SITE)  =  10  •  P_code(l) + P_code(2) 
j_CODE(i_SrrE)  =  P_index(  PARM_code(LSITE)  +  1 ) 

END  DO 

IF  (  j_CODE(l)  .NE.  j_CODE(2)  )  THEN 
WRnE(6,*) 'ERR;  The  PARMAMETER  values  differ  between' 
WRITE(6,*)  ’  tiie  two  ionosonde  input  files!' 

WRITE(6,*) 'ACT:  Execution  has  been  halted.' 


A-22 


WRrTE(6,*)  'Please  check  these  files  and  try  again.' 

STOP 

ELSE 

LCODE-j_CODE(l) 

IF 

'WRITE(6,*)  •MONTH; i_MONTH.’  PARM:  ^CODE 

If  not  all  of  the  parameters  have  yet  been  read  from  the  current  file 
AND  the  current  parameter  code  exceeds  the  code  counter  "N^CODE",  then 
set  "N^CODE”  equal  to  the  parameter  CODE  counter  "i_CODE''.  This  state 
implies  that  not  all  of  a  month's  data  blocks  have  yet  been  input  and 
processing  of  the  monthly  data  cannot  proceed. 

Otherwise,  determine  if  all  of  the  data  blocks  necessary  to  process 
the  current  month's  data  have  been  read  for  both  stations.  If  so, 
process  the  month's  data. 

IF  ( i_MONTH  .EQ.  MONTH_l  )  THEN 

IF  ( i_CODE  .LE.  N_CODE  )  THEN 

"LCODE  <  N^CODE"  means  that  N__CODE  has  been  reached  a  maximum  from 
a  first  pass  through  all  blocks  in  a  previous  month's-worth  of  data. 

In  this  case,  the  first  month's  data  has  been  read  and  the  first  block 
in  the  second  month's  blocks  has  been  detected.  Process  the  first  months 
data 

DOLDAY=l,31 
DO  i_HOUR«  1, 24 

CALL  HFLINK2(  i_DAY,  i_HOUR,  LINK  flag,  - 
MONTHLYJlag) 

END  DO 
END  DO 

IF  ( month_OLD  .2Q.  0  )  monthjOLD  =  i_MONTH  - 1 
IF  (  month_OLD  .EQ.  0  )  monthjOLD  =  12 

IF  (  MONTHL Y_flag  )  THEN 

CALL  NRMLYZE(  i_YEAR,  month_OLD,  MONTHLYJlag  ) 
END  IF 

Njnonths  =  N  months  +  1 
month_OLD  ^ImONTH 
FIRST  month  .FALSE. 

OLD jCfONTH  =  MONTH(l) 


ELSE 


N_CODE»LCODE 
END  IF 

ELSE  IF  ( i_MONTH  NE.  MONTO,!  )  THEN 

DOLDAY-1,31 
DOi  HOUR- 1,24 

CALL  HFLINK2(  i  DAY.  i  HOUR.  LINK  flag,  - 
MONTHLY  flag) 

END  DO 
END  DO 

IF  (  mondi.OLD  .EQ.  0  )  month_OLD  -  i_MONTH  - 1 
IF  (  month^OLD  .EQ.  0  )  month_OLD  -  12 

IF  (MONTHLY  flag)THEN 

CALL  NRMLYZE(  i  YEAR,  month.OLD,  MONTHLY_flag  ) 
ENDIF 

N^montfis  —  Njnonths  + 1 
month.OLD  « i_MONTH 
F^T  ffloath  -  JALSE. 

OLDJMONTH  -  MONTH(l) 

ENDIF 

"  If  die  cunrent  parameter  code  is  not  to  be  used,  read  the  next  data 
"  block  and  return  to  start  of  loop. 

CODE^flag  -  .TRUE. 

IF  ( .NOT.  CODE  maskC  i  CODE  )  )  THEN 
DOi-1,31 

ERR_LINE(1)  «  ERR_LINE(1)  +  1 
READ(4,702,END-200)  RECORDJ(l) 

ERR_LINE(2)  -  ERR_LINE(2)  +  1 
READ(7,702,END-200)  RECORDJ(2) 

END  DO 

CODE  flag -.FALSE. 

ENDIF 

IF  ( .NOT.  CODE Jlag  )  CYCLE 

"  Determine  the  upper  and  lower  bin  limits  of  the  parameter  for  either 
"  frequency  (iJDIM  -  1)  or  height  (i_piM  -  2) 
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i  DIM-P  type(i  CODE) 

NJ)iiis  -  TO^ins^LDIM  ) 

Compute  parameter  increment  value 

d  P-(P  W  CODE)-PJow(i  CODE))- 
/  FLOAT(  ifPbinsCLCODE) ) 

START  loop  for  each  day  of  the  current  mon& 

i_DAY«0 

DAY.flag-.TRUE. 

DO  WHILE  ( i_DAY  U.  31  .AND.  DAYJlag  ) 
LDAY-LDAY+1 

Read  next  input  record  and  divide  into  24  5>character  strings 
Check  the  first  20  characters  to  make  sure  th^  are  nonblank.  If  they 
are  blank,  then  a  new  record  should  be  read. 

ERR  LINE(1)  •  ERR  LINE(1)  +  1 
READ(4,702,END-W)  RECORD  1(1) 

ERRJLINE(2)  »  ERR  LINE(2)  + 1 
READC7,702,END-997)  RECORD J(2) 

WRITE(6,*)  ’SriE  1  INPUT  RECORD’ 

WRITE(6,*)  RECORD  1(1) 

wRrrE(6,*)  'sriB  2  input  record* 

WRITE(6,*)  RECORD J(2) 

Check  to  SM  if  current  day  is  widiin  interval  to  be  processed. 

IF  (  ( i  YEAR  .EQ.  YEAR  i  .AND.  i  MONTH  .EQ.  MONTHJ  - 
.Al^.  i  DAY  .LT.  DAY  i )  .OR  - 
.NOT.  DAY_mask(i J)AY)  )  CYCLE 

IF  ( i  YEAR  .EQ.  YEAR  f  .AND.  i  MONTH  .EQ.  MONTH  f - 
.A^D.  i  DAY  .GT.  DAYJ ) 

DAYJlag  “.FALSE. 

READ  flag  “.FALSE. 

CYCLE 
END  IF 

START  loop  for  each  hour  of  the  current  day. 

HOUR_flag*.TRUE. 


iHOUR-0 

DO  WHE.E  ( i_HOUR  XT.  24  .AND.  HOUR_flag  ) 
i_HOUR-LHOUR+l 
N_dy_hr(LHOUR)  -  N_dy_hr(i_HOUR)  +  1.0 

NJir*N3r+  1-0 

"  Qieck  to  see  if  current  hour  is  within  interval  to  be  processed 


IF  ( ( i  YEAR  .EQ.  YEARJ  .AND.  i  MONTH  EQ.  MONTHJ  - 
.AND.  i  DAY  XQ.  DAY_i  .AND.  iJHOUR  .LT.  HOURJ )  - 
.OR  .NOT.  HOUR  mask(i  HOUR)  )  CYCLE 


"  IF  ( i  YEAR  .EQ.  YEAR  f  .AND.  LMONTH  EQ.  MONTH  f  - 

"  .Al^.  i  DAY  XQ.  DAY  f  .AND.  i  HOUR  .GT.  HOURJ )  THEN 

"  HOUR  flag -.FALSE. 

"  DAYjiag- JALSE. 

*  Rept  flag  -  .FALSE. 

"  CYOX 

•  END  IF 

"  Determine  die  parameter  value  for  each  ionosonde  station. 

k_srrE»o 
DOi_SITE-l,2 
P_value  =  0.0 

k_CHAR- 5  *  ( iJIOUR-  1 )  +  1 

IF  (RECORD  I(i  SITE)(k  CHARk_CHAR44) - 
.EQ.'  ’)CYCXE 

"  Read  scaled  parameter  for  the  current  hour  at  the  ionosonde. 

P  char(l)  =  RECORD  I(LSrrE)(k_CHARk_CHAR) 

Plchar(2)  =  RECORD J(LSITE)(k_CHAR+l:k_CHAR+l) 

P  char(3)- RECORD J0_SrrE)(k_CHAR+2:k_^CHAR+2) 

CLchar  -  RECORD  J(T_SrrE)(k_CHAR+3  ;k_CHAR+3) 

D_char  -  RECORD  J(LSrrE)^_CHAR+4;k_CHAR+4) 

IF  (P_char(l)  .NE. "  )  THEN 

P  l-REAL(ICHAR(P_char(l))-48) 

P” 2»REAL(ICHAR(P_char(2))-48) 


P J  -  REAL(  ICHAR(  P_char(3)  )  -  48 ) 

P~flag  -  .TRUE. 

ELSE  IF  (  P_char(2)  .NE.  " )  THEN 
P  1-0.0 

P“ 2-REAL(ICHAR(P  char(2))-48) 

PJ  -  REAL(  ICHAR(  Plchar(3)  )  -  48  ) 
P_flag».TRUE. 

ELSE  IF  (  P_char(3)  NE. "  )  THEN 

PJ-0.0 

PI2-O.O 

P  3  -  REAL(  ICHAR(  P_char(3)  )  -  48  ) 
P_flag-.TRUE. 

ELSE 

P  value -0.0 
PJlag-. FALSE. 

END  IF 

IF(P_flag)THEN 

P  value  «  (  100.0*P  1  +  10.0*P  2  +  P_3  )  - 
/P_scale(iJ>^ 

k_SITE»k_SITE  +  1 

CALL  P ARMBN2(  i_SnE,  i_DAY,  i_HOUR,  i_CODE,  - 
P_value,  Q_char,  D jchar  ) 

END  IF 
End  of  station  loop. 

JF  (  k_SITE  .EQ.  2  )  N_GOOD  =  N_GOOD  +  1.0 
END  DO 

End  of  hourly  loop  for  current  day,  parameter,  month,  and  year  (file). 
END  DO 


End  of  duly  loop  for  current  parameter,  month,  and  yeu  (file). 


END  DO 


”  End  of  current  parameter  &.  month  loop  for  year  (file). 

END  DO 

*  All  input  records  in  die  current  input  file  (or  year)  have  been  read. 

200  CLOS'  iT»4) 

”  Input  results  from  next  year  (or  file) 

END  DO 

"  Compute  HF  links  for  last  month 

DOi  DAY- 1,31 
DOi  HOUR- 1,24 

CALL  HFLINK2(  i  DAY.  i_HOUR,  LINKJlag,  MONTHLY_flag  ) 
ENDDO 
END  DO 

Njmondis  -  N_months  +  1 
IF  (MONTHLY  flag)  THEN 

CALL  NRh£YZE(  LYEAR,  month^OLD,  MONTHLY^flag  ) 

END  IF 

mondi  OLD“i_MONTH 
FIRSTimonth- JALSE. 

”  Ou^nit  die  probability  distribudons  for  each  parameter  code,  both 
"  for  each  Imur  and  ovts  all  hours  in  a  column  format 

IF  ( .NOT.  MONTHLYJlag  )  THEN 

"  All  input  files  have  been  read.  Normalize  each  distribution  count 
"  by  the  numbers  of  samples  for  each  paramater  code. 

CALL  NRMLYZE(  i^YEAR,  LMONTH,  MONTHLY_flag  ) 

WRrrE(3,704)  STATION  iiame(l),  STATION_name(2) 

704  FORMAT(/10x,TONOSO^E  PARAMETER  STATISTICS  AT ',  - 

/10x,a20,*&',a20) 

WRnE(3,705)  HOUR_i,  DAY_i,  MONTHJ,  YEARJ,  HOUR_f,  DAY_f,- 
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MONTHJ.  YEARJ 

705  FORMAT(/10x,'IONOSONDE  measurements  processed  for  the  period;  - 

/lOx,  Trom:  3(i2;.  to  \3(i2.\  •).i4) 

"  Output  HF  link  propagation  characteristics. 

WRITE(3,730) 

730  FORMAT(//10x,'HF  LINK  PROPAGATION  VERSUS  HOUR  &  2-HOP  MODE') 
"  Output  hourly  results  for  each  two-hop  propagation  mode 
DOi_MODE=l,4 
DOj_MODE“1.4 

WRITE(3,734) -PROPAGATION  MODE: 

MODE_name(i_MODEJ_MODE).  - 
(i.i-1,24) 

WRITE(3,736)  '100  Pr*.  - 

(100.0*N  mode  hr  dy(i  MODEj_MODE,i_HOUR),i  HOUR*l,24),  - 
lOO.O^N  mc^e  dy(i_MODEJ_MODE) 

WRITE(3,737)  'AVGFRQ',  - 

(AVG  fMHz  mode  hr(i  MODEJ  MODE,i  HOUR),  - 
i  HOUR«l,24),  AVG  fMHz  mod<i  MODEJ  MODE) 

WRITE(3,737)  ’SDVFRQ',  - 

(SDV  £MHz_modeJir(LMODEJ_MODE,i_HOUR),  - 
i_HOUR=l,24),  SDV  fMHz  mode(LMODEJ_MODE) 

"  \VRITE(3,737)'AVGANG’,- 

"  (AVG  adg_mode  hr(i  MODEJ_MODE,i_HOUR),  - 

"  i  HOUR=l,24),  AVG  adg_mode(i_MODEj  MODE) 

"  WRrrE(3,737)  'AVGDST, 

"  (AVG_dkm  modc_hr(i  MODEJ_MODE,i_HOUR),  - 

"  iJHOUR=l,24),  AVG_dkm_mode(i_MODEJ_MODE) 

END  DO 

END  DO 

'WRrrE(3,732) 

WRITE(3,738)  2.0*Ds_kin, '  FREQ '(MHZ) ',  ( i,  i  *  1, 24  ) 

DOi_DISP=l,N_DISP 

WRITE(3,742)  f  RANGE(i  DISP),  - 

(  100.0*fhfflz_hr_disSLHOUR,i_DISP),  i_HOUR=  1 , 24  ),  - 
100*fh5lz_disp(i_DISP) 


END  DO 


WRrnE(6,*)  •"GOOD"  hour  percentage.  N_GOOD  /  N_hr 
WWTE(6.*)  Tercentage  of  NO  LOS  hours;  N^noLOS  /  N_hr 

WRITER, 746) 

WRITE(3,748)  2.0*Ds_km.  'ANGLE ',  ( i,  i  -  1. 24  ) 
DOi_ANG-O.N_ANG 

a__dg  ■  adg^low  +  i__ANG  •  d_a 

WRITE(3.752)  a_dg,(adg_hr_bin(iJiOlIR,LANG),LHOUR=l,24),  - 
adg_bm(i_ANG) 


END  DO 

WWTE(3.754)  (adg_hr  bin(i  HOUR,N_ANG^-l),LHOUR»1.24),  - 
adgJ)m(NJ^G+r) 

WRITE(3,756) 

WRITE(3.758)  2.0*Ds_kin. '  RANGE',  ( i.  i  =  1. 24  ) 


DOi_DST-0,N_DST 
dkm  =  dkm Jow  +  i_HGT  *  d.dkm 
WRITE(3,762)dkni,- 

(dkmJir_binCi_HOUR,i_DST),i_HOUR»l,24),. 

dkmJbinO.DST) 

ENDDO 

WRITE(3,764)  (dkm  hr  bin(i_HOUR,N_DST+l),i_HOUR=l,24),  - 
dlan_bm(NjbST+l) 


End  of  distance  loop 
END  IF 

CLOSE(UNIT=2) 

IF  (  MONTHLY_flag  )  THEN 
CLOSE(UNIT=8) 

py  jgp 

cxosE(UNrr»3) 

END  IF 

STOP 
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997  WRITE(6,*)  'Error  on  input,  probably  wrong  number  of  - 
characters  in  rea>rd' 

CLOSE(UNIT«2) 

CLOSE(UNrr-3) 

STOP 

99i  WR]TE(6,*)  'Error  on  input,  probably  no  access  to  input  file' 
CLOSE(UNrr«2) 

CLOSE(UNrr-3) 

STOP 

999  WRITE(6,*)  'No  ionosonde  input  files  were  listed' 

CLOSE(UNIT-2) 

CLOSE(UNIT-3) 

STOP 

600  FORMAT(a20) 

601  FORMAT(a25,lx.a25) 

700  FORMAT(a20,lx,a20) 

701  FORMAT(lx,'i[nput  filename:  ',a30) 

702  FORMAT(al20) 

734  FORMAT(//lx,  al9,  aS,  33x,  llour  of  the  day  - 
/7x,24(lx,i3,  lx), 'NET') 

736  FORMAT(  lx,  a6, 25(f4.1,lx)  ) 

737  FORMAT(lx,a6,25(f4.1,lx)) 

732  FORMAT(//10x,'*****  MAXIMUM  USABLE  FREQUENCY  **♦**’) 

738  FORMAT(//5x,'  Link  Range;  ',f7.1,'KM'- 

//Ix,  ii6, 48x,'H0ur  of  die  day  (UT)'  > 

/Ix,a5,24(2x,i3),'  NET') 

742  FORMAT(  lx.  a6, 25(f4.1,lx)  ) 

746  FORMATfl'/ZlOx,'*****  RAY  ELEVATION  ANGLE  ****♦') 

748  FORMAT(//5x, '  Link  Range:  '.H.l.'KM'- 

//Ix,  96, 48x,  'Hour  of  the  day*  - 
/7x,24(i3,2x),'NET') 

752  FORMAT(  lx,  f6.2, 25f5.3  ) 

754  FORMAT(  lx, '  >  ',  25f5.3  ) 

756  FORMATCl'/ZlOx,'*****  FIRST-HOP  PATH  LENGTH  **♦**') 

764  FORMAT(lx,'  >  '.  25f5.3  ) 

758  FORMAT(//5x, '  Link  Range;  '.fZ.l.'KM*- 

//Ix,  a6, 48x,  'Hour  of  the  day*  - 
/7x,24(i3,2x),'NET') 

762  FORMAT(lx,f6.1,25£5.3) 


END 


SUBROUTINE  NRMLYZE(  LYEAR,  LMONTH,  MONTHLY_flag  ) 

*  Declare  amys  for  HF  link  calculations 

REALM  IMHzJow,  ACHzJigh,  dkmjow,  dkmjigh 
REALM  AiHzJu_disp,  fKfflz_di^ 

REAL*4  ad^kfw,  adgjb^  d_fMHz,  d_dkm 

REALM  ajimit,  adgj>in,  adgju’_bin,  adg_modeJbin 

REALM  adgjnodejhr^bin,  AVG_adg_mode_hr,  AVG_adg_mode 

REALM  delimit,  dkmjbin,  dkni_lir_bin,  dkmjnodeJ)in 

REALM  dkmjnodejirjbin,  AVG_dkm_inode_lir,  AVG_dkm_mode 

REALM  ^limh,  fMHz_bin,  £MHz_lir_bin,  fMHzjnode_bin 

REALM  fMHzjnodeJhr  bin,  AVG  fMHz  modejir,  AVG  fMHzjmode 

REAL*4  SDv3fHzjnode_hr,  SDV^fMIfajnode,  N_dy_hr,  N_hr,  N_GOOD 

REAL*4  N_dy,  NJ»r_dy,  N_inodc__dy,  N_mode_hr_dy,  N_smp,  N_noLOS 

1NTE^R*4  i  DAY.  i  HOUR,  i  DST,  i  FRQ,  i  MODE,  i  SITE 
INTEGERM  iT.DST,  NJFRQ,  n'aNG  ” 

INTECXRM  NJbins,  month.OLD 

CHARACTER*SL  lieading(3),  MODE  name(l;4,l:4) 

C31ARACTER*25  DIST.heading 

"  Declare  variables  f(»rMUFcalcuiadons 

REALM  in_MUF,b>IUF,DUM  1,DUM  2,N_km  hr,N  km 
REALM  MUF J0(1 :24),  MUF_50(1 :24),  MUF_9b(l-24) 

REALM  ANGJ0(1:24),  ANG_50(1:24),  ANG_90(1:24) 

INTE(®RM  i_YEAR,  i_MONTH.  N_DISP 

LOGICAL*!  P10_FLAG,  P50_FLAG,  P90J1AG,  MONTHLYJlag,  NONZERO_flag 
LOGICAL*!  TESTJFLAG 
CHARACTER*25  LNK_file 
”  D^ine  all  common  blodcs. 

COMMON/IONSNP/foF2(2,31.24),  M3000F2(2,31,24),  hF2(2,31,24),  - 
foFl(2,31,24),  M3000F1(2,3I,24),  hFI(2,31,24),  - 
hF(2,31,24).foE(2,31,24),- 
hE(2,31,24).  foEs(2,31,24),  fbEs(2.31,24),  - 
fEs(2,31,24),liEs(2,3I,24) 

"  finin(2,3 1,24),  fxl(2,3 1,24) 
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COMMON  /MISC/  Rejcm,  Rearth_km 


COMMON  yUMIT/  fMHz_low,  fMHz_hgh,  N__FRQ,  dkm^low,  (lkin_hgh,  N  DST,  - 
adg_low,  N_ANG,  LNK_file,  d_fMHz,  - 

d.adg,  As_rd,  Ds_km 

•  COMMON  /ANGLE/  adg_bin(0;61).  - 

•  adg_lir_bin(l:24.0:61X  «d8^mode_bin(l:4, 1:4,0:61),  - 

"  adgjnodeJ»rJ»in(l  :4,1 :4,1 :24,0:61),  - 

"  AVG__adg_mile_hr(l:4,l:4,l:24),  AVG  adg_mode(l:4,l:4), - 

"  ANG_hr_frq(l:24,0:61,l:2) 

COMMON  /ANGLE/  ANG_hr_frq(l  :24,0:61.1:2) 

"  COMMON  /RANGE/  dkm_biii(0:61),  - 

dkin^hrjbm(l:24,0:61),  dkm_mode_bm(l:4, 1:4,0:61),  - 
dkm_modeJhr  bm(l  :4, 1:4, 1:24,0:61),- 
"  AVG_dkm3n^c_hr(l:4,l:4.1:24),  AVG_dkm_mode(l:4,l:4) 

COMMON  /FRQBN/  f_liimt(l:12,l:2),  fMH2_bm(0:61),  - 

fMH2_hrJbin(l:24,0:61),  fMHzjnodc  bm(l:4,l:4,0:61),  - 
fMH2jnodcJ»rJ)m(l  :4,1 :4,1 :24,0:61)7- 
AVGjfMHzjBiodcJir(l:4, 1:4,1:24),  - 
SDV_fMH2lnodc3r(l:4,l:4,l:24),  - 

AVGJMH2jnode(l:4,l:4),  SDV  fMHz  mode(l:4,l:4),- 
fMHzJir_disp(l  .24,1:12),  fMHzjdisp(l:r2) 

COMMON /COUNT/ N_dy,N  hr  dy(l:24),N_mode  dy(l:4,l  :4),N  GOOD,- 
Njnode_hr_dy(l:4,l^,l;24),  N_dyjii(l:24),  NJir,  - 
N^noLOS 

DATAN_piSP/12/ 

IF  (N_dy  XT.  0.5  )  RETURN 

DOi  FRQ«0,N  FRQfl 

fll!H2_biii(i  ^Q)  «  £MHz_bm(LFRQ)  /  N  dy 

ENDDO 

"  DOi_ANG=l,N  ANG 

•  ad£.bin(i  ANG)”®  adgj)m(i  ANG)  /  N_dy 

•  ENDDO 

“  DOi  DST=1,N  DST 

-  dkm  bm(i  DST^  *  dim  biii(i_DST)  /  N  dy 

•  ENDDO 


D0LDISP-1,N  DISP 
fNfflz  disp(i  DISP)«fMH2_disp(i  DISP)/N  hr 
END  DO 

DOtHOUR-1,24 
N_8mp  -  N_dy_hr(LHOUR) 

IF  (  N_smp  .LT.  0.5  )  CYCLE 
DOiJDISP-l,N_DISP 

iMHz  hr  disp(i  HOUR,i  DISP)«fMHz  hr  disp(i  HOlJR,i  DISP)/N  smp 
ENDDO  --- 

N_smp  »  N_hr_dy(i_HOUR) 

IF  (  N_smp  XT.  0.5  )  CYCLE 

DO  i.FRQ  -  0,  N_FRQ+1 

aiH2_hrJ)mOJIODA,i  FRQ)-fMH2  hr  bm(i_HOUR,i_FRQ)  /  N  smp 
IF  ( iJFRQ  .GT.  0  .AND.  MONTHLY  )  THEN 
£MHz.hrJ)m(i  HOUR,*  FRQ)»£hffl2  hr  bin(i  HOUR,i  FRQ)- 
+  fMHz  hr  bin(i  HOUR,r  FRQ^l) 

ENDIF 

ENDDO 


DOLANG»l,N  ANG 

adgJjr_biii(i_HdUR,i_ANG)  =  adg_hr_bin(tHOUR,i_ANG)  - 
/N  smp 

ENDDO 


II 

It 


DOi_pST=l,N  DST 

dlcm_hr_biii(LHOUR,i_DST)  «  dkm_hr J)mCi_HOUR.i_DST)  - 

/Njanp 


ENDDO 


DOi_MODE*=l,4 

DOj_MODE*l,4 

"  Compute  the  PDFs  for  each  bin,  mode,  and  hour 

N  8mp  =  N  mode  hr  dy(i_MODEJ_MODE»i_HOUR.) 

BF” [  N_smp  XT.  07s  fCYCLE 

N  mode  hr  dy(i_MODEJ_MODE,LHOUR)  = - 

Njnode_hr_dy(i_MODEJ_MODE,i_HOUR)  /  N_dy 
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D0LFRQ-1,N  FRQ 

fMH2^mode_,hr_bm(i_MODEj^MODE,i_HOUR,LFRQ) « - 
fMH2_mode_^_biii(i_MbDEJ_MODE,i_HOUR,i_FRQ)  - 
/N  smp 

ENDDO 

"  DOLANG*l.N  ANG 

•  «d£inode_br  bin(i  MODEJ  MODE,i  HOUR,LANG)«- 

"  tdg^mode_br_bin(i_MODEj^MODE,i_HOUR,LANG)  - 

•  /N  smp 

"  ENDDO 

"  DOLDST- l.N  DST 

"  dkm_modc_br_bin(i  MODEJ  MODE.i_HOUR,i_DST)  « - 

“  <ikmjnodc_brbiii(i_MODEJ_MODE,i_HOUR,i_DST)  - 

"  /N  smp 

"  ENDDO 

"  Compute  tbe  statistics  for  each  mode 

SUM  fMHz«AVG  fMHz  mode  br(i  MODEJ  MODE,i  HOUR) 
SD^ JMHz  -  SDVJ>dH2jn<5eJff(i_MODij_MODE,i_HOUR) 

AVG^fMHz  «  SUM  fMHz  /  N  smp 

AVG JMHzjnode JirCi JMODE J_M0DE,LH0UR)  «  AVG JMHz 

IF  (  N_smp  .GT.  1.5 )  THEN 

SDVJMHz  -  N  smp  ♦  SDV2  £MHz  -  SUMJMHz**2 
IF  (  SDV  fMHzXT.  0.0  )  SDV  fMHz  -  0.0 
SDVJbfflz  -  SQRT(  SDy_fMHz  /  (  N  smp  *  (  N_smp  - 1.0  )  )  ) 
SDV_fMHzjnodeJir(iJM6DEJ_MOra,i_HOUR)  *  SDV^fMHz 
ENDIF 

SUM_adg  «  AVG_adgjaode_br(i_MODEJ_MODE,i_HOUR) 

AVG  adg^SUM  adg/N  smp 

AVGladg_mode_hr(i_MODEJ_MODE,iJIOUR)  «  AVG.adg 

SUM_dkm  =  AVG_dkm_mode_br(LMODEJ_MODE,i_HOUR) 

'  AVG  dkm  =  SUM_dkm  /N  smp 

•  AVGldkmjnodeJirCi_MODEJ_MODE,i_HOUR)  =  AVG_dkm 
ENDDO 

ENDDO 


ENDDO 


DOLMODE-1.4 

DOLMODE-1,4 

N  smp-N  mode_dy(i  MODEj_MODE) 

IF  (  N_smp  .LT.  0.5  )  CYCLE 

N  mode  dy(i  MODEJ  MODE)-- 

Njmode_dy(LMODEj^MODE)  /  N Jy 

SUM  fMHz-AVG  fMH2_mode(i_MODEJ_MODE) 

SDV2_fMH2  -  SDV_fMHz_mode(LMODEj_MODE) 

AVG  fMHz-SUM  fMHz/N  smp 
AVGlfMHz_modc(rMODEJ  JdODE)  -  AVG.fMHz 

IF(N  smp.GT.  1.5)THEN 
SDV  fMHz-  n  smp  *  SDV2  fMHz  -  SUM  fMHz**2 
IF  (  SDV  fMHz'iT.  0.0  )  SDV  fMHz  -  0.0 
SDV  fMHz- SQRT( SDV  fMHz/(N  smp ‘(N  smp- 1.0))) 
SDV" fMHz  mode(i  MOdIj  MODE) -SDV  fN^ 

ENDIF 

SUM  adg-AVG  ad&.mode(i  MODEj_MODE) 

AVG  adg-SUM  adg/N  smp 
AVGladgLmode(i_MOD^j_MODE)  «  AVG_adg 

SUM  dkm-AVG  dkm  mode(i_MODEJ.MODE) 

AVG" dkm  -  SUM  dkm7  N  smp 
AVGldkm_modcCi>fODEJJ40DE)  -  AVGjIkm 

DOi  FRQ»1,N  FRQ 

fWHz_mode_bin(LMODE  j  J40DE,i JFRQ)  -  - 

£MHzjniodcJbiii(i_MODEJ_MODE,i_FRQ)  /  Njanp 

ENDDO 

DOi  ANG»1,N  ANG 
ad£modeJ>in(i_MODEJ_MODE,i_ANG)  -  - 

adg  mode  bin(i_  MODEJ_MODE»i_ANG)  /  N_smp 

ENDDO 

DOi_pST-l,N_DST 
dkm_modeJbm(i_MODEJ_MODE,i_pST)  =  - 

rflcm  modeJbiii(i_MODEJ_MODE,i_DST) / N_smp 

ENDDO 
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END  DO 


END  DO 

Exit  if  statistics  are  not  being  collected  monthly. 

Output  percentage  of  "good”  hours  in  which  a  usable  path  was  found. 

IF  ( .NOT.  MONTHLY^flag )  RETURN 

Compute  tenth,  fiftieth,  and  ninetieth  percentiles  for  the  link  MUF 
add  results  to  the  LNK^file 

DOi.HOUR*  1,24 

TEST  FLAG -.FALSE. 

IF  ( i_HOUR  .EQ.  10  .AND.  i_MONTH  .EQ.  2  )  TEST^FLAG  =  TRUE. 

DO  i_FRQ  -  0,  N_FRQ+1 

DUM  2»fMH2  hr  binO  HOUR.i  FRQ) 

DUM_1  -  fMH2lhrIbin(rH0UR,NJFRQ+l) 

IF  (DUM  1  .GT.  0.0)  THEN 
fMHz  hr  bin(i  HOUR,!  FRQ)  =  DUM  2  /  DUM  1 
ELSE 

fMHz  hr_bin(i  HOUR,i  FRQ) -0.0 
ENDIF 

f » fMHzJow  +  FLOAT( i_FRQ-l)  *  dJMHz  +  0.5  ♦  dJMHz 
END  DO 

PIO  FLAG -.FALSE. 

P50  FLAG -.FALSE. 

P^_FLAG-. FALSE. 

DOi_FRQ=l,N_FRQ-l 

x2  =  fMH2jirJbin(i_HOUR,iJFRQ+l) 

IF  (  x2  .GT.  0. 10  .AND.  .NOT.  P10_FLAG  )  THEN 

P10_FLAG  =  .TRUE. 

xl  -fMHz.hr  bin(i  HOUR,LFRQ) 
yl  -  fMHzJow  +  (  FLOAT(  i_FRQ  -  1 )  +  0.5  )  *  dJMHz 


y2  ■  yl  +  d_fMHz 
DUMJ«x2-xl 

IF(DUM  1  .NE.O.O)THEN 
m  MUF-(y2-yl)/DUM  1 
b  MUF-(yl  •x2-y2*xl)/DUM  1 
0UM  2-in  MUF*0.10  +  b_MUF 
ELSE  ~ 

DUM  2»fMHz  low  +  (FLOAT(i  FRQ- 1  )  +  0.5 ) *  d  fMHz 
END  IF 

IF  ( DUM  2  .GE.  fMHzJow )  THEN 
MUF  10(i  HOUR)* DUM  2 
ELSE 

MUF  10(i  HOUR)  *  fMHz  low  +  0.5  •  d.fMHz 
ENDff 

DUM_1  -  ANG_hr_frq(i_HOUR,LFRQ.l) 

IF  (  DUM_1  .GT.  0.0  )  THEN 

DUM^2-ANG  hr  frq(i  HOUR,!  FRQ,2) 

ANGJO(i_HOim)  *  DIM_2  /  DUMJ 

ELSE 

"  Detomine  average  angle  below  current  MUFIO  frequency 

DUMJ  *0.0 
DUM  2*0.0 

NON^O  flag*. FALSE. 

DOj_FRQ«i  FRQ-1, 1,-1 
IF(NON2mO  flag) THEN 
CYCLE 
ELSE 

DUMJ  *  ANGJirJrq(iJIOURj_FRQ,l) 

IF  (  DUM_1  .ME.  0.0  )  THEN 
NONZ^O  flag  *  .TRUE. 

DUMJ  *  ANG_hr_frq(LHOUP*J_FRQ.2) 

END  if” 

END  IF 
END  DO 

"  Determine  average  angle  above  current  MUFIO  frequency 

DUM  3  *  0.0 
DUMJ  *0.0 
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NONZERO  flag  -  .FALSE. 

DOj  FRQ-i  FRQ+l.N  FRQ 
IF(NONZ^O  flag)'raEN 
CYCLE 
ELSE 

DUM  3  -  ANG  hr  frq(i  HOURJ  FRQ.l) 

IF  (DUM  3  .NE.  0.0 ) TOEN 
NONZl^O  flag -  .TRUE. 

DUM  4 -ANG  hr  frq(i  HOURJ  FRQ,2) 

endif” 

END  IF 
END  DO 

DUM  5 -DUM  1+DUM  3 
IF(DUM  5.NR0.0)THH^ 

DUM  6-DUM  2  +  DUM  4 
ANG  10(i  HOUR)»DUNr6/DUM  5 
ENDIF 

END  IF 

ENDIF 

IF  (  x2  .GT.  0,50  .AND.  .NOT.  P50_FLAG  )  THEN 

P50JFLAG-.TRUE. 

xl  -iMHz  hr  biii(i  HOUR,i  FRQ) 
yl  -  fMHzJow  +  (  FLOAT(  FfRQ  -  1 )  +  0.5  )  *  d^fMHz 
y2  -  yl  +  d_fMHz 

DUM_l-x2-xl 

IF  (DUM  1  .NE.  0.0)  THEN 
m  MUF  =  (y2-yl)/DUM  1 
b  MUF-(yl*x2-y2*xl)/DUMJ 
DUM_2  =  m_MUF  *  0.50  +  b_MUF 

■PT  RP 

DUM_2  =  fMHzJow  +  (  FLOAT(  i_FRQ  -  1 )  +  0.5  )  ♦  dJMHz 
ENDIF 

IF  (DUM  2  XSE.  fMHzJow  )  THEN 
MUF  50(i  HOUR)»DUM_2 
ELSE  *" 

MUF_50(i  HOUR)  =  fMHzJow  +  0.5  *  d_fMHz 
ENDIF 


DUM.l  -  ANG.hr Jrq(i_HOUR,tFRQ,l) 
IF(DUM_1  .GT.  0.0)  THEN 


DUM_2-ANG.hr.frq(i  HOUR,i.FRQ,2) 
ANG_50(i_HOim)  -  DljM.2  /  DUM.l 

ELSE 

Determine  avenge  angle  below  current  MUFSO  frequency 

DUM.l -0.0 

DUmI2»0.0 

NON^O  flag  -  .FALSE. 

DOj  FRQ-i  FRQ-l,  1,-1 
IF  TnONZ^O  flag  )  THEN 
CYCLE 
ELSE 

DUM  1  -ANG  hr  fiq(i.HOURj_FRQ,l) 
IF  (DUM  1  .ME.  o7o )  THEN 
NONZ^O  flag  -  .TRUE. 

DUM  2 -ANG  hr  frq(i  HOURJ  FRQ,2) 
END  IF 
END  IF 
END  DO 

Detennine  average  angle  above  current  MUFSO  frequency 

DUM  3-0.0 
DUM.4-0.0 

NON^O  flag  -  .FALSE. 

DOj  FRQ-i.FRQ+l,N  FRQ 
IF  CNONZERO  flag  )  TOEN 
CYCLE 
ELSE 

DUM.3-ANG  hr.frqO.HOURj  FRQ,1) 
IF  (DUM  3  .NE.  0.0 ) THEN 
NONZ^O.flag  =  TRUE. 

DUM  4 -ANG  hr.frq(i.HOURJ  FRQ,2) 
END  IF 
END  IF 
END  DO 

DUM  5  -  DUM.l  +  DUM_3 
IF(DUM  5.NR0.0)TH^ 

DUM  6"^  DUM  2  + DUM  4 
ANG”50(i  HOUR) -DUM  6 /DUM  5 
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END  IF 


END  IF 
END  IF 

IF  (  x2  GT.  0.90  .AND.  .NOT.  P90 )  THEN 

P90_FLAG- .TRUE. 

xl  -fMHz_hr_biii(i  HOUR,i_FRQ) 
yl  -  fMHzJow  +1 FLOATX  i_FRQ  -  1 )  +  0.5 )  •  d_fMHz 
y2  •  yl  +  d_IMHz 

DUMJ-X2-X1 

IF  (  DUM_1  .NE.  0.0  )  THEN 
m_MUF  *  (  y2  -  yl  )  /  DUM  1 
b3fUF-(yl  ♦x2-y2*xl')/DUM  1 
DUM  2»m_MUF*0.90  +  b  MUF  ” 

ELSE 

DUM  2«£MHz  low  +  (FLOAT(i  FRQ- 1  )  +  0.5 ) •  d  fMHz 
ENDIF 

IF  (DUM  2. GE.  fMHz  low)THEN 
MUF  90(i  HOUR)«DUM  2 
ELSE 

MUF_90(^HOUR)  »  fMHzJow  +  0.5  ♦  dJMHz 
ENDIF 

DUM_1 «  ANGJirJrq(LHOUR,i_FRQ,l) 

IF  (  DUM_1  .GT.  0.0  )  THEN 

DUM  2  =  ANG  hr  fi:q(LHOUR,i  FRQ,2) 

ANG_90(i_HOi:m)  «  DUM_2  /  DUM J 

ELSE 

Detennine  average  angle  below  current  MUF90  frequency 

DUMJ«0.0 

DUM2-0.0 

NON^OJlag  =  .FALSE. 

DOj  FRQ-i  FRQ-1, 1,-1 
IF  ^NONZ^O  flag  )  THEN 
CYCLE 


ELSE 

DUM  1  -  ANG  hr  frq(i  H0URJ_FRQ,1) 
ff(0UM  l.NE.O.O)TOEN 
NONzAo  flag -.true. 

DUM  2->WG  hr  frq(i  HOURJ  FRQ,2) 

END  IF 
END  IF 
END  DO 

*  Determine  average  angle  above  cuirent  MUF90  frequency 

DUM_3-0.0 

DUMjl-0.0 

NON^O  flag -.FALSE. 

DOLFRQ-i  FRQ+1,N  FRQ 
IfTnONZ^O  flag)THEN 
CYCLE 
ELSE 

DUM.3  -  ANG  hr  frq(i  HOURj  FRQ,1) 

IF  (DUM  3  .NE.  0.0)*raEN 
NONZ^O  flag -.TRUE. 

DUM  4-^G  hr  firq(i  HOURj  FRQ,2) 

END  IF 
END  IF 
END  DO 

DUM  5 -DUM  1+DUM  3 
IF  (  DUM  J  .NeT  0.0  )  THH^ 

DUM  6-DUM  2  +  DUM  4 
ANG”90(i  HOUR)  =  DUm16/DUM_5 
ENDIF 

ENDff 

ENDIF 

IF  (  P10_FLAG  AND.  P50_FLAG  .AND.  P90_FLAG  )  CYCLE 
END  DO 
END  DO 

"  Output  results  by  ^)pending  to  disk  file 
"  OUTPUT  MUF  values  exceeded  90%  of  the  time 

IF  ( i  MONTH  .EQ.  0  )  i  MONTH  -  MONTH_l 
WRITE(8,700)  i_YEAR,  i_MONTH,  ( i,  i  -  0. 24  ) 
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700  FORMAT(/2x.  12x, '  10%  NOT  EXCEEDED  (FOT)',  lx,  i4,  lx,  i2,  - 
/2x,  6x,  2S(  lx,  i3,  lx ) ) 

WRITE(8,702)  TOT:  MUF  10(24),  - 

(MUF  10(i  HOUR),i  HOUR=l,24),  - 
•ANG:  ',ANG  10(24),  - 
(ANG  10(i  HOUR),i_HOUR»l,24) 

702  FORMAT(  lx,  a6, 25f5.1  /lx,  a6, 25f5.1 ) 

”  OUTPUT  MUF  values  exceeded  50%  of  the  time 

WRITE(8,704)(i,i-0.24) 

704  FORMAT(2x.  12x. '  50%  NOT  EXCEEDED  (MUF)',  - 
/2x,6x.25(  lx,i3,  lx)) 

WRITE(8,702)  -MUF;  ',  MUF  50(24),  - 

(MUF  50(i  HOUR),i*HOUR»1.24).  - 
•ANG:  •.  ANG  50(24),  - 
(ANG_50(i_HOUR),iJ!OUR»l,24) 

*  OUTPUT  MUF  values  exceeded  10%  of  the  time 

W!UTE(8,708)(i,i-0,24) 

708  FORMAT(2x.  12x, '  90%  NOT  EXCEEDED  (HPF)',  - 
/2x,6x,25(  lx,  i3,  lx)) 

WWTE(8,702)'HPF:  '.MUF  90(24),- 

(MUF  90(i  HOUR),rHOUR-l,24),- 
'/WG:  •,ANG  90(24).- 
(AI  _90(i  HOUR),!  HOUR=l,24) 


adgj>m  =  0.0 
adgjirjbin  =  0.0 
adgjmodejbin  =  0.0 
adgjnode_hrJbin  “  0.0 
AVG_adgjnodeJir  =  0.0 
AVGjidgjnode  =  0.0 
ANG_hrJRrq  »  0.0 
dkmjbin  «  0.0 
rflrm  hr  bin  ^  0.0 
dkm_modeJbin  “  0.0 
d]cm_mode_hr_bin  =  0.0 
AVG_dkm_mode_hr  *  0.0 
AVGjlkmjmode  **  0.0 
j|^4Hz_bm  *  0.0 
a4Hz_hr_bin  *  0.0 
AiHz_mode_bin  =  0.0 
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IMHzjDiode_hr_bin  ■  0.0 
AVG  fMHz  hr  »  0.0 

SDVJMHzjnode Jit  -  0.0 
AVGJMHz  mode* 0.0 
SDVJMHzjnode -0.0 
fMHz.hr.disp  -  0.0 
iMHzJisp  *  0.0 
N_dy-0.0 
NJhr jly  -  0.0 
Njnode jly  «  0.0 
Njnode Jn_dy  -  0.0 

MUFJO-0.0 

MUFJO-0.0 

mufIoo-o.o 

ANGJO-0.0 

ANG_50-0.0 

ANG_90-0.0 

RETURN 

END 


"  Subroutine  BEARING 

II 

"  Purpose:  Compute  the  latitude  and  longitude  of  two  locations 

*  in  radians,  tire  bearing  angles  from  true  north  and 

"  the  Great  Circle  path  distance. 

SUBROUTINE  BEARING(  STATlON^name,  COLAT JXj,  LNG J>G,  Ds_km,  - 
As jrd,  azimutiijdg,  MONTHLY Jlag  ) 

R£AL*4  LAT  dg(2),  LNG  dg(2),  Dsjon,  As  dg,  a^uth  dg(2),  - 
S  LAT(2)7c  LAT(2),  T_LAT(2),  S  LNG(2),  C_LNG(2), - 
T^ LNG(2),  B  dg(2),  FLAG(2),  S_dLNG(2),  C^dLNG,  DUM,  - 
cbLAT_dg(2),  Re_km.  Rearth_to 

CHARACTER*20  STATION_name(2) 

LOGICAL*!  MONTHLY  Jlag 

COMMON  /MISC/  Re_km,  Rearth^km 

*  Define  constants 

pi-4.0*ATAN(1.0) 
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nldg>  180.0 /pi 
dgrd  •  pi  / 180.0 


Loop  on  «ich  ionosonde  station  location 
DOi_SrrE-l,2 

LAT_dg(LSITE)  -  90.0  -  COLAT_dg(tSITE) 

S_LAT(i  SITE)  -  SIN(  dgrd  •  LAT  dg(i  SITE)  ) 

C.LAT(i  SITE)  -  COS(  dgrd  *  LAT  dg(i  SITE)  ) 

T J,AT(i_SnE)  -  TAN(  dgrd  •  LAT_dg(i_SnE)  ) 

S  LNG(i  SITE)  -  SIN(  dgrd  ♦  LNG  dg(i  SITE)  ) 

C_LNG(i  SITE) «  COS(  dgrd  •  LNG  dg(i_SITE) ) 
T_LNG(i_SnE)  «  TAN(  dgrd  •  LNG_dg(i_SITE)  ) 

END  DO 

C  dLNG  «  COS(  dgrd  ♦  (  LNG  dg(l)  -  LNG  dg(2)  )  ) 

S  dLNG(l)  -  SIN(  dgrd  *  (  LNG  dg(2)  -  LNG  dg(l)  )  ) 

S.dLNG(2)  -  SIN(  dgrd  •  (  LNGJg(l)  -  LNGldg(2)  )  ) 

Path  length  in  Idlometers  and  Great  Circle  arc  length 

DUM  «  C_LAT(1)  *  C_LAT(2)  •  C_dLNG + S_LAT(1)  •  S_LAT(2) 
Ds  Jon  *  Rearthjcm  *  ACOS(  DUM  ) 

As_rd  »  Ds  Jem  /  Rearth Jem 

• 

Bearing  angles  fiom  1  to  2  and  2  to  1  measured  clockwise  from 
true  north 


DOLSITE=I,2 
j  SITE** 2 

W  (  LSHE  .EQ.  2  )LSITE  =  I 

DUM  =  C  LAT(i  SITE)  *  T.LATO.SITE)  -  S_LAT(LSnE)  *  C_dLNG 
B_dg(i_SnE)  =  rddg  *  ATAN2(  S_dLNG(tSITE),  DIM  ) 

FLAGCLSriE)  =  SIGN(  1.0,  S_dLNG(tSnE) )  *  SIGN(  1.0,  B_dg(i_SnE) ) 

IF  (  FLAG(i_SITE)  .GT.  0.0  )  THEN 
azimuth  dg(i  SITE)  =  B_d^i__SITE) 

ELSE 

azimuth  dg(i  SITE)  =  B_dg(i_SnE)  +  180.0 
ENDF 


IF  (  «imuth_dg(tSITE)  XT.  0.0  )  THEN 
aziinuth  d^i  SITE)  *  azimuth  dg(i  SITE) +  360.0 
ENDIF  “  ' 

END  DO 

WRITE(6.700) 

WIUTE(6,701) 

WRITE(6,703)  ( i,  STATION  iiame(i).  COLAT_dg(i).  LNGJg(i).  i  =  1. 2  ) 
WRITE(6,705)  (  azimuth  dg(i).  i  - 1, 2  ).  - 
2.0 *Ds  km.2.0*Ds  km/ 1.609 
WRITE(6,700) 

IF  (  MONTHL Y_£lag  )  THEN 

WRnE(8,700) 

WWTE(8,701) 

WRrrE(8,703)  ( i.  STATION_iiame(i),  COLAT  dg(i).  LNG  dg(iX  - 
i-1.2) 

WRITE(8,705)  (  azimuth  dg(i).  i  *  1, 2  ),  - 
2.0  *  Ds  km,  2,0*  Ds  km  / 1.609 
WRITE(8,700)  ’ 

ELSE 

WRITE(3,700) 

WRITE(3,701) 

WRrrE(3,703)  ( i,  STATION_name(i).  COLAT  dg(i),  LNG_dg(i),  - 
i-1,2) 

WRrTE(3,705)  (  azimuth^dgO),  i  - 1, 2  ),  - 
2.0  *  Ds  km,  2.0  *  Ds^km  / 1.609 
WWTE(3,700)  ” 

ENDIF 

700  FORMATC//*— - 

r - •) 

701  FORMATCy  lONOSONDE  STATION  GEOGRAPHY ') 

703  FORMAT(y  Site i2, 3x,  a20 ,  - 

/I  lx,'  CoLatitude:  ',f6.1,'  DEG - 
/I  lx,'  East  Longitude:  'i'S.  1/  DEG ') 

705  FORMAT(y  Site  l->Site2  Site2->Site  1',- 

y  Bearing  angles:  ',f6.1,'DEG  ',5x,f6.1,'DEG’,- 
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t  Link  Great  Circle  path  length:  ',f7.2,'  KM*.  - 
t  or’,n.2,'MI*) 

RETURN 

END 


AJUt  HFLINK2.F0R  Source  Program 


SFREEFORM 

SLARGE 


*  Subroutine  HFLINK2 

*  Purpose:  Process  lONOSONDE  data  files  from  two  stations  following  the 

*  conventions  established  in  UAG-23  and  UAG-23A  and  detennine 

*  the  optimum  frequencies  for  all  possible  two-hop  HF  radio  links 

”  with  reflection  points  above  each  ionosonde. 

SUBROUTINE  HFLINK2(  i_DAY.  i_HOUR,  LINK.flag,  MONTHLY_flag  ) 

"  Declare  variable  types  and  dimension  arrays 

"  Declare  arrays  for  HF  link  calculations 

REAL*4  m>fUF.  b.MUF,  DUM^l,  DUM_2,  NJmJir,  N.km 

REALM  fMHzJow,  fMHzJhgh,  dkmjow,  dlmjbgh 

REALM  fMHzJujdisp,  fMHzjiisp 

REALM  adg^low,  adgjhgh,  d_fMHz,  d_dkm 

REALM  ajimit,  adgj)in,  adgjuj>in,  adgjnodejbin 

REALM  ac^m^ejhr.bin,  AVGjsdgjnodeJir,  AVGjsdgjnode 

REALM  djimit,  dkmjbin,  dkmjirjbin,  dkmjnodejbin 

REALM  dkm^modejirjbin,  AVG_dkm_mode_hr,  AVG_dkm_mode 

REALM  fjimit,  fMHzJbin,  fMHz_hrJbin,  £MHz_modeJbin 

REALM  fMHz_modeJirJbin,  AVG_fMHzjnodeJir,  AVG_fMHzjnode 

REALM  SDV_fMHz_mode_hr,  SDVJMHzjnode,  NjaoLOS 

REAL*4  N__dy,  N_hr_dy,  Njmodc_dy,  N_modeJu_dy,  N_dyjir,  NJir,  N_GOOD 

REALM  LAT_dg(2),  LNG_dg(2),  Dsjon,  azimudi_dg(2X  Rmirthjcm 

REALM  hv  lm(4,4^Xfo  MHz(4,4,2),  Re  km^alfa  rd(2),MUF  MHz(2) 

REALM  f  MUF  MHz(4,4  d_MUF_km(4,4,2),  a_MUF_dg(4.4) 

REAL  h  UNK  ^  f_LINK_MHz,  a_LINK  dg,  D  km,  ANGJirJirq 
REAL*4'mUF  fctr,  hF2_km,  hFl.km,  hE  1^  h£s  km 
REALM  Kf;  k'VALUE,  r  km(2),  phi  rd(2),  TO  rd(2) 

REALM  foF2,  ficF2,  fiF2.  M3000F2,  hF2,  hp£2,  foFl,  fxFl,  - 
M3000F1,  hFl,  hF,  foE,  foE2,  hE,  hE2,  foEs,  - 
fxEs,  fbEs,  fEs,  hEs,  foFldS,  finin,  M3000Fld5,  hfldS,  - 
foL  fid,  fini,  12000,  n,  Lcxx 

EOTGERM  i_DAY,  iJiOUR,  i J)ST,  i_FRQ,  i J40DE,  i_SITE 

INTECffiRM  N_DST,  N_FRQ,  N_ANG 

INTEGERMN  bins,  month_OLD 

INTEGERM  i  YEAR,  i  MONTH 

INTEGERM  ONE,  NJWSP 
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CHARACTER*5L  heading(3),  MODE  naine(l;4.1;4) 
CHARACTER*25  LNK  file 
CHARACTER*25  DIST_hcading 

LOGICAL*!  PIO  FLAG,  P50  FLAG.  P90  FLAG,  MONTHLY  flag 
LOGICAL*!  Fa£  flag,  MODE  flag(4,4.2).  NON2XRO  flag 
LOGICAL*!  INTV,  URN,  D  flag,  LINK_flag(4.4).  GOOD  LINK 
LOGICAL*!  TEST  FLAG 


Define  common  blocks 

COMMON/IONSNP/foF2(2.3!,24),  M3000F2(2,3!.24).  liF2(2,31,24),  - 
foF!(2,3!,24).  M3000F  1(2,3!, 24),  hF!(2,31,24),  - 
hF(2.31.24).foE(2,31.24).. 

IiE(2.31,24).  foEs(2.31,24).  fl)Es(2.31,24).  - 
fEs(2,31,24),  hEs(2,31.24) 

finin(2,31,24).  1x1(2,31,24) 

COMMON  /MISC/  Re_km,  Rearth^km 

COMMON  /LIMIT/  fMHzJow,  fMHzJigh,  N_FRQ,  dkmjlow,  dkmjigh,  N_DST,  - 
adgjow,  adgji^  N_ANG,  LNK^file,  djfMHz,  d^dbn,  - 
d_adg,  As_rd,  Ds Jon 

COMMON  /ANGLE/  adgj»in(0:61),  - 

adgjurjbin(l:24,0;61),  ad^mode_bin(l:4,!;4,0:6!),  - 
adg_modeju_bin(l:4,!:4,!;24,0:61),  - 
AVG_adg:  mode  hi(  1 :4, !  :4, 1  ;24),  AVG_adgjnode(l  :4, 1 :4),  - 
ANG_hr Jrq(!  ;24,0:6!,  1:2) 

COMMON  /ANGLE/  ANG_hr_frq(!;24,0:6!,!;2) 

COMMON  /RANGE/  dkm_bin(0;6!),  - 

dkm_hrjbin(!:24,0:6!),  dlm_modejbin(!:4,l;4,0;61),  - 
dkm  mode_hr  bin(l:4,l:4,l;24,0:6!), - 
AVG_dkm_modeJir(l:4,l:4,l;24),AVG_dkm_mode(l;4,l;4) 

COMMON /FRQBN/f_limit(l:  12, 1:2),  fMHzJ)in(0:61),  - 

£MHzJjrJ)in(l:24,0:61),  £MHzjnode_bm(l:4, 1:4,0:61),  - 
A(iHz_mode  hr  bin(l:4,l:4,l:24,0:61), - 
AVG  fMHzIm^e_hr(l:4,l:4,l:24),- 
SDV  IMHzjnode  hr(l:4,l:4,l:24). - 

AVG  fMH2_mode(l:4,l:4),  SDV_fMHz_mode(l:4,l:4),  - 
fMHThr  disp(l:24,l:12).fMH2_disp(l:12) 


COMMON  /COUNT/  N_dy.  N_hr_dy(l  ;24),  N_mode_dy(l  :4. 1 :4),  N_GOOD.  - 
N  mode  hr  dy(  1:4, 1:4, 1:24),  N  dy_hr(l:24),  N_hr,  - 
N  noLOS  “ 


DATAN  DISP/12/ 


"  Compute  constants 

S_3000  -  SIN(  3000.0  /  (  2.0  •  Rearth_km ) ) 

C_3000  -  COS(  3000.0  /  (  2.0  *  Rearth_km ) ) 
pi-4.0*  ATAN(l.O) 
rddg- 180.0 /pi 
dgrd  -  pi  / 180.0 

hv_km  =  -1.0 
fo_MH2--1.0 
f  MUF_MHz--1.0 
aJLINlCdg^O.O 
a  MUF  dg-0.0 
dlLINlfkm-O.O 
d_MUF_km  -  0.0 

MODEJlag-. FALSE. 

LINK Jlag-. FALSE. 

"  DETERMINE  REFLECTION  HEIGHTS  FOR  EACH  AVAILABLE  MODE 
DOi_HOP-l,2 

"  WRITE(6,*)  i  YEAR,  i_MONTH,  iJDAY,  iJiOUR,  i  HOP 
"  WRITE(6,*)  •foF2: ',  foF2(i_HOP,i_DAY,i_HOUR) 

IF  (  foF2(LHOP,i_DAY.i_HOUR)  .GT.  0.0  )  THEN 

ff  (  M3000F2(i  HOP,i  DAY,i_HOUR)  GT.  0.0  )  THEN 
MUF  fctr«1^000F2(i_HOP,LDAY,i_HOUR) 

M3060F2(i  HOP,i_DAY,i_HOUR)  =  -1.0 
DUM  1  -MUF  fctr/k  VALUE( 3000.0 ) 

DUM”  1  -  SQRf(DUM_l**2  -  1 ) 
hF2  to  -  Rearth_to  ♦  J S_3000  /  DUM_1  +  C_3000  -1.0) 
ELSEIf  (  hF2(i  HOP,i  DAY,i_HOUR)  GT.  0.0  )  THEN 
hF2_to  -  hF2(i_HbP,i_DAY,i_HOUR) 
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hF2(i  HOP.i_DAY,i_HOUR)--1.0 
ELSE 

hF2  km  »  325.0 
ENDff 

DOLMODE«l,4 

IF(i  HOP.EQ.  1)THEN 

MODE  nag(lj_MODE.l)-.TRUE. 
fo  J_MODE,l)  «  foF2(l,i_DAY,i_HOUR) 
hv  km(lj_MODE.l)-hF2_km 
END  IF 

IF(i  HOP.EQ  2)'raEN 

”MODE_flag(j_MODE,1.2)  » .TRUE. 
fo_MHzO_MODE,l,2)  -  foF2(2,i_DAY,LHOUR) 
hv  kmO  MODE,l,2)®  hF2_kin 

endif“ 

END  DO 

foF2(i_HOP,iJ)AY,LHOUR)  *  -1.0 
ELSE 

DOj_MODE-l,4 
IF  ( i  HOP  .EQ.  1 )  THEN 

“mode  flag(lj_MODE,l)». false. 
fo_MH2(ljlMODE,l)  »  -1.0 
hv_km(lJ_MODE,l)  *  -1.0 
ENDIF 

IF(i_HOP.EQ.2)THEN 

MODEJlagO_MODE,l,2)=  JALSE. 
fo  MHzO_MODE,l,2)*-1.0 
hvJanO_MODE,l,2)  »  -1.0 
ENDff 

END  DO 

ENDIF 

WRrTE(6,*) '  foFl : foFl(i_HOP,i_DAY,i_HOUR) 

IF  (  foFl(i_HOP,i_DAY,i_HOUR)  .GT.  0.0  )  THEN 


IF(M3000Fl(i  DAY.i  HOUR)  .GT.  0.0 )  THEN 
MUF  fctr«I^000FT(i  HOP.i  DAY.i  HOUR) 
M300bFl(i  HOP.i  DAY.i  HOUR)«.r.O 
DUM  1-^  f^/k  VALUE( 3000.0) 

DUM  1  *  SQRf(  DUM  1**2  - 1 ) 

DUM_2  -  S JOOO  /  DUMJ  +  C JOOO  -  1.0 
hFl  km*Reartb  km* (DUM  2) 

ELSEIF ( UFIO  HOP,i_DAY.i  HOUR)  GT. 0.0 ) THEN 
hFl  km-tfl(i  HOP,i  DAY.i  HOUR) 
hFl(i  HOP,i  DAY.i  HOUR)* -1.0 
ELSE  IF  ( hF(i  HOP.i  DAY.i  HOUR)  GT.  0.0  )  THEN 
hFl^km  -  hF(i  HOP.LDAY.LHOUR) 
hF(i  HOP.i  DAY.i  HOUR)  =  -1.0 
ELSE 

hFl_km  «  200.0 
ENDff 

DOj_MODE=1.4 

IF(tHOP  JEQ.  1)THEN 

MODE_flag(2j_MODE.l)  =  .TRUE, 
fo  MHz(2J  MODE,l)  =  foFl(l.i  DAY,i_HOUR) 
hvJcmCZ  j_MODE.  1)  =  hFl_km 
END  IF 

IF(i  HOP.EQ.2)THEN 

”MODE_flagO  MODE,2,2)  =  .TRUE, 
fo_MHz(j_MODE,2,2)  =  foFl(2,i_DAY,i_HOUR) 
hv  km(j  MODE,2,2)  =  hFl_km 

endif” 

END  DO 

foFl(i_HOP,LDAY,i_HOUR)  =  -1.0 
ELSE 

DOj_MODE=l,4 

IF(i_HOP.EQ.  1)1HEN 

MODE_flag(2j_MODE,l)  =  .FALSE, 
fo  MHz(2j  MODE,1)  =  -1.0 
hv_km(2j_MODE,l)  =  -1.0 
ENDff 

IF(i  HOP.EQ.2)THEN 

MODE_flag(j_MODE.2.2)  =  .FALSE. 


fo^MH2(j^MODE.2,2)  -  -1 .0 
hv  kin(j  MODE,2,2)*-1.0 
ENDIF 

END  DO 

ENDIF 

WRITE(6,*)  •  foE:  *,  foE(i_HOP,LDAy.i_HOUR) 

IF  ( foE(i_HOP,LDAY,LHOUR)  GT.  0.0  )  THEN 

IF  ( hE(i  HOP,i  DAY.i  HOUR)  .GT.  0.0  )  THEN 
hE  lcin«hE(i  HOP.i  DAY.i  HOUR) 
hE(i  HOP,i  DAY,i  HOUR)- -1.0 
ELSE 

hE  km  -  110.0 

endIt 

DOLMODE-1,4 

IF(i  HOP.EQ.  1)THEN 

MODE  flag(3j  MODE,l)«  TRUE. 
fo_MHz(3j  MODE,l)-foE(l,i  DAY,!  HOUR) 
hv  km(3J  MODE,l)-hE  km 
ENDIF 

ff  (i_HOPm2)THEN 

MODE_flagO_MODE,3,2)  -  .TRUE, 
fo  MHzO_MODE,3,2)-foE(2,i  DAY,i_HOUR) 
hv_km(j_MODE,3.2)  =  hE^tai 
ENDIF 

END  DO 

foE(i_HOP,iJ)AY,i_HOUR)  -  -1.0 
ELSE 

DOj_MODE-l,4 

IF(i  HOP  .EQ.  1 ) THEN 

"mode  flag(3J_MODE,l)-. FALSE, 
fo  MHz(3  J_M0DE,  1)  -  -1 .0 
hv_km(3  J_MODE,  1)  -  -1 .0 
ENDIF 


IF(i  H0P.EQ.2)THEN 

MODE  flag(j  MODE,3,2)«  FALSE. 
fo,MH2(j_MODE3.2)  «  -1.0 
hv  km(j  MODE,3,2)--1.0 
END  IF 

END  DO 

END  IF 

WRrrE(6.*)  •  foEs;  •.  foEs(i_HOP.i ' 

IF  ( foE5(LHOP,iJDAY.i_HOUR)  .GT.  0.0 )  THEN 

IF  (  hEs(iJHOP,i J>AY,i  HOUR)  .GT.  0.0  )  THEN 
hEs.fcBi»hEs(i  HOP,rDAY,i  HOUR) 
hEs(i  HOP,i  DAY,i  HOUR)*”  1.0 
ELSE 

hEs  km  *  120.0 
ENDff 

DOj_MODE»l,4 

ff  (LHOPEQ.  1)THEN 

MODE_flag(4j  MODE,l)*  .TRUE. 
foJv!Hz(4J_MdDE,l)  *  foEs(l,i  DAY,i_HOUR) 
hv_km(4J_MODE,l)  «  hEsJm 
END  IF 

IF(i_HOP.EQ.2)THEN 

MODE Jlag(j>fODE.4,2)  =  .TRUE, 
foj^filzb  MdDE,4,2)«foEs(2,iJDAY,i  HOUR) 
hv  km(j_MODE,4,2)*hEs_km 
ENDF  ” 

END  DO 

foEs(LHOP,i_DAY,i_HOUR)  *  -1.0 
ELSE 

DOj_MODE*l,4 

IF(i  HOP.EQ.  1)THEN 

MODE  flag(4J_MODE,l)  =  .FALSE. 
fo_MHz(4jlMODE,l)  *  -1.0 
hv  km(4j  MODE,1)«-1.0 


END  IF 


ff  (LH0P.EQ.2)THEN 

MODE  flagO  MODE.4.2)*’ .FALSE. 
fo.MHrO  MODEA2)  « -1.0 

hv_kmO_MODE,4,2)  -  -1 .0 
END  IF 

END  DO 

END  IF 

"  WRITE(6.*) 

"  WRITE(6,777)  'fo  MH2(l.*,i  HOP):  ’.(fo  HOP),i=l,4) 

"  WRITE(6.777)  To  MHz(2,*,i_HOP):  ‘.(fo  MH2(2,i.i_HOP),i*l,4) 

*  WRITE(6.777)  To^MHz(3,*,i_HOP):  •,(fo_MHz(3.i,i_HOP),i=l,4) 

“  WRrrE(6,777)To  MHz(4,*.i  HOP);  '.(fo  MH2(4,i,i_HOP),i«l,4) 

"  WRITE(6,*) 

"  WWTE(6,777)  •hv_,km(l,*,i_HOP);  '.(hv  HOP),i*l,4) 

"  WRITE(6,777)  Tiv  kiii(2.*.i  HOP):  ’.(hv  km(2,i,i  HOP),i*l,4) 

*  WWTE(6,777)  •hvjan(3,*,i  HOP):  '.(hv  km(3,i,i  HOP),i=l,4) 

"  WRrrE(6,777)  liv_km(4,*,i  HOP):  \(hv  km(4,i,i  HOP),i*l,4) 

"  WWTE(6,777) 

"  WjaTE(6,*)  •MODE_flag(l,*,i_HOP):  ’.(MODE  flag(l,i,i  HOP),i=1.4) 

“  WRITE(6,*)  •MODE_flag(2,*,i  HOP):  ‘.(MODE  flag(2,i,i  HOP),i«l,4) 

"  WRITE(6,*)  'MODE_flag(3,*,i  HOP):  ’.(MODE  flag(3,i,i  HOP),i»l,4) 

"  WWTE(6,*)'MODE_flag(4,*,i  HOP):’,(MODE  flag(4,U_HOP),i=l,4) 

"  WRITE(6,*) 

ENDDO 

777  FORMAT(lx,a20,4(f6. 1, lx)) 

”  Exit  if  no  two-hop  propagation  is  possible. 

K  COUNTED 
DOi  MODE*  1,4 
DOj  MODE  =1,4 

LH^  flag(i  MODEJ_MODE)  =  MODE_flag(i  MODEJ_MODE,l)- 
”  .AND.  MODE_flag(i_MODEj_MODE,2) 

IF  ( LINK_flag(i  MODEJ_MODE) )  K_COUNT  =  K_COUNT  +  1 
ENDDO 
ENDDO 

"  IF  (  K_COUNT  .GT.  0.0  )  N_GOOD  =  N_GOOD  +  1 .0 
-  WRITE(6,*)  LINK_flag(l,*):  ',(LINK.flag(l,i),i=l,4) 
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"  WWTE(6,*)  ^INK^fl«g(2.•):  ‘.(LINK  flag(2,i),i»1.4) 

"  WRITE(6,*)  'LINK_flag(3,*):  ’.(LINK  nag(3,i).i«l,4) 

"  WIUTE(6,*)  •LINK_fl^(4,*);  •,(LINK_flag(4,i).j«1.4) 

*  WRrrE(6,*) 

*  Detomine  the  propagation  range  achieved  for  each  reflection  height. 

f_LINK_MHz»  10000.0 
MUF_hto--1.0 
GOODJLINK  « .FALSE. 

N,LOS«0 

IX)i_MODE*1.4 

DOj_MODE-l,4 

Drop  out  of  loop  if  current  two-hop  mode  does  not  have  adequate  data. 
IF  ( .NOT.  LINKJlag(LMODEj_MODE)  )  CYCLE 
N_LOS»NJLOS+l 

“  Determine  if  overl)^)ping  radio  LOS  is  adueved  ^ven  the  two  virtual 

*  reflection  hdghts. 

rjan(l)  -  Rearth_km  +  hvJan(i_MODEj_MODE,l) 
aI_LOS_rd  *=  ACOS(  Reaithjkm  /  r_km(l)  ) 

rjm(2)  »  Rearth_km  +  hv_km(i_MODE  j_MODE,2) 
a2_LOS jrd  =  ACOS(  Rearthjcm  /  r Jan(2)  ) 

IF  (  al_LOS_rd  +  a2_LOS  rd  iT.  As  rd  )  THEN 
LINK_flag(i  MODEJ_MODE)  =  .FALSE. 

N  LOS*N  LOS-1 
ENDff 

IF(  .NOT.  LINKJlag(LMODEj_MODE)  )  CYCLE 

"  Detemune  the  ground  reflection  point  where  the  angle  of  incidence  is 
"  equal  to  tiie  angle  of  reflection.  This  calculation  requires  root  finding 
"  and  employs  a  modified  Newton's  method. 

a^LOWj-d  =  0.0 

IF  (  a2JLOS_rd  .LT.  As_rd  )  a  LOW  rd  =  As  rd  -  a2_LOS_rd 
a_HGH_rd  =  aLLOS  rd 
IF  (  al_LOS_rd  .GT.  As_rd  )  a_HGH_rd  =  As_rd 
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theu  rd  -  ROOT(  t  LOW  rd,  a_HGH  rd.  r  km,  As  rd,  INTV,  ITRN ) 
IF  ( .NOT.  INTV  OR.  NOT.  ITRN  )  THEN 
LINK_flag(LMODEj_MODE)=  FALSE. 

END  IF 

IF  ( .NOT.  LINK_flag(i_MODE  j_MODE) )  CYCLE 

alfa_rd(l)-2.0*thcta_id 
alfa_rd(2)  *  2.0  •  As_rd  -  alfa_rd(l) 


GOOD_LINK».TRUE. 

DO  i_HOP  -1.2 

DUM_1  -  SIN(  alfa_rd(LHOP)  /  2.0  ) 

DUM  2  -  1  +  hv  km(i  MODEJ_MODE,i  HOP)  /  Rearth  km  - 
-  COS(  aIfa_rd(CHOP)  /  2.0  ) 

IF  (  DUM_2  IT.  1.0  )  THEN 

IF  (  1.0e38  *  DUM  2  .GT.  DUM  1  )THEN 
phi  rd(i  HOP)  -  ATAN2(  DUM  I,  DUM  2) 

ELSE 

phi  rd(i  HOP)- pi/ 2.0 

endIf 

ELSE 

phi_rd(i_HOP)  -  ATAN2(  DUM_1,  DUM_2  ) 

ENDIF 

TO_rd(i_HOP)  =  ( pi  -  alfa_rd(iJIOP) )  /  2.0  -  phi_rd(i_HOP) 


Compute  the  "secant  corrected"  factor  k. 


d_GND  km  =  Rearth  km  ♦  alfa_rd(i_HOP) 

IF  ( iJHOP  .EQ.  1  fd  MUF  km(i_MODEj_MODE,l)  -  d_GND_km 
Kf-  k_yALUE(  d_GND_km  ) 


Compute  the  maximum  usable  frequency  for  the  the  current  hop. 

MUF  MHz(i  HOP)  =  Kf*fo  MH2(i_MODEj_MODE,i_HOP) - 
“  /COS(phi_rd(i_HOP)) 


END  DO 
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IF  (MUF_MH2(1)  XT.  0.0  .OR.  MUF_MH2(2)  XT.  0.0 )  CYCLE 

*  Two>hop  MUF  is  the  mimmum  of  the  MUFs  of  the  two  hops  and  the 
"  two-hop  TO  angle  is  given  by  the  corresponding  value  of  alfa. 

"  WRrrE(6,*)  *As  dg;  rddg*As  rd,'  theta  dg: rddg*theta  rd 
"  WRrrE(6,*)  'alfa  dg(l): rddg*alfa  rd(l).  - 

*  •  alfa  dg(2); rddg*alfa  rd(2r 

"  WRITECd/rphi  dg(l):',rddg*phi_rd(l).. 

"  'phi  dg(2);',rddg*phi  rd(2) 

"  WRrrE(6,*)  TO  dg(l): rddg*TO  rd(l),  - 

"  ’TO  dg(2):',rddg*TO  rd(2) 

"  WRrrE(6.*)‘ - • 

IF(MUF  MHz(l) .LT. MUF  MH2(2))THEN 
f  MUFjMHzCi  MODEJ  MODE) -MUF  MH2(1) 
a  MUF  dg(i  MODEJ  MODE)  »  rddg  •  fo_rd(l) 

ELSE 

fMUF  MHz(i  MODEJ  MODE) -MUF  MHz(2) 
a  MUF  dg(i  MODEJ  MODE)  -  rddg  *  fo_rd(2) 

E^^IF 

”  The  overall  MUF  is  the  minimum  of  all  the  usable  MUF  values. 

IF  ( f  MUF  MHz(i  MODEJ  MODE)  .LT.  f  LINK  MHz  )  THEN 
f  Liric_MIfa-f_MUF_MHz(i  MODEJ  MODE) 
a_LINK  dg-a  MUF  d^i  MCn)EJ_MC^E) 
d  LINK“km-d  MUF  km(i  MODEJ_MODE,l) 

EbTOIF 

END  DO 

END  DO 

IF  (  fLINK_MHz  .GT.  1000.0  )  GOOD_LINK  -  .FALSE. 

"  WRrrE(6,*) 

"  WRITE(6,*)  'GOOD  LINK.  GOOD.LINK 

"  WIUTE(6,777)  'f  MUF_MHz(l  J):  (f_MUF_MHz(l J)J=1,4) 

"  WRITE(6,777)  •f_MUFJ^2J): (f_MUF_MHz(2J)J-l,4) 

"  WRrrE(6,777)  'f_MUF_MHz(3J): (f_MUF_MHz(3J)J-l,4) 

"  WRITE(6,777)  ’rMUF>lHz(4J): (f_MUF_MHz(4J)J-l,4) 

"  WRITE(6,*) 'fLINK_MHz:  ',fJLINK_MHz 

"  INCREMENT  THE  APPROPRIATE  NUMBER  DENSITY  BINS. 

"  Exit  if  no  link  was  possible. 
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IF  ( .NOT  GOOD  LINK  )  THEN 
IF  (  N  LOS  .LE.  0  )  N_noLOS  »  N_noLOS  +  1 .0 
f  LWK  MHz-O.o" 

^TUKN 
END  IF 

Total  number  of  hours  with  usable  data,  the  frequencies,  angles,  and 
distanMs. 

N.GOOD  -  N_GOOD  +  1.0 
Njdy  “N^dy  + 1.0 

CALL  BINS  0(  f  LINK  MHz,  fMHzJow,  fMHz.hgh,  fMHzJin,  N_FRQ  ) 
CALL  BINSl0(  a.LINi^dg,  adg_.low.  ad&_hgh,  adg_bin,  nJaNG  ) 

CALL  BINS_0(  d_LINKJan,  dkmjow,  dlcm Jigh,  dkmj>in,  N_DST  ) 

Increment  the  number  of  days  for  which  data  is  av^able  in  the  current 
hour. 

N_hr_dy(i_HOUR)  -  NJir_dy(i_HOUR)  + 1.0 

CALL  BINS  1(  i  HOUR,  f  LINK  MHz,  fMHz  low,  fMHz  hgh,  24,  - 
fMHzJirJbin,N_FRQ) 

CALL  BINS  1(  i_HOUR,  a_LINK_dg,  adgjow,  adg_hgh,  24,  - 
adgJir_bin,N  ANG) 

CALL  BINS_1(  i_HOUR,  dLLINK^km,  dkmjow,  dkmjigh,  24,  - 
dkm JirJ)in,  NJ)ST  ) 

Incremoit  the  number  of  days  for  \^ch  data  is  available  in  the  current 
hour  for  each  available  two-hop  propagation  mode. 

DOi_MODE=l,4 

D0LM0DE=1,4 

IF  ( .NOT.  LINK Jlag(i_MODE  j_MODE)  )  CYCLE 

N_mode_dy(iJMODEj_MODE)  *  N_mode_dy(i_MODEj_MODE)  +  1.0 
N_modeJlir  dy(i  MODEJ_MODE,iJIOUR)®- 

Njnode_hr_dy(iJvlODEj_MODE,LHOUR)  + 1.0 

AVG_fMHz_mode_hr(i_MODE jJ40DE,i JIOUR)  *  - 
AVG_fMHz_mo«ie_hr(i_MODEJ_MODE,i_HOUR)  +  f_LINK_MHz 


DUM-f  LINK  MHz* f  LINK  MHz 
SDV.fKfflz.mo^  hr(i  MODEJ  MODE,i  HOUR)«- 

®VJMHzjn^eJr(LMODEj  MODE,LHOUR)  +  DUM 


"  IF  ( i  MODE  £Q.  1  .AND.  j  MODE  .EQ.  1  AND.  i_HOUR  EQ.  1 )  THEN 
"  WRnE(6,*)  i  MODE,  j  MODE,  i  HOUR,  f  LINK  MHz 

•  WRITE(6,*)  'AVG.  ’,  AVG  fMHz  mode  hr(i  MODE j_MODE,i  HOUR) 

•  WRnE(6,*)  'SDV: ',  SDV  fMHz  mode  hr(i  MODEj_MODE.i  HOUR) 

•  ENDIF 

•  AVG  zdg_mode  Mi  MODEJ  MODE,!  HOUR)-- 

"  AVG_Kla_m^e  hr(i  MODEJ  MODE,!  HOUR)  +  a  LINK  dg 

"  AVG.dkm  mode  Mi  MODEJ  MODE,!  HOUR)- - 
“  AVG_im_m^ejMLMODEJ_MODE,i_HOUR)  +  dJLINK^km 

AVGJMHz  mode(i  MODEJ  MODE)-- 

AVG  fMHz_mode(i  MODEJ  MODE)  +  f  LINK  MHz 
SDV.fMHz  mode(i  MODEj~  MODE)  - - 

SDV  fMHz  mode(i  MODEJ  MODE)  +  DUM 

•  AVG_adg_m^eCi  MODEJ“mODE)  -  - 

•  AVG  adg^mode(i  MODEJ_MODE)  a  LINK  dg 

"  AVG  dkm  m^e(i  MOD^  MODE)-- 

"  AVG_5an_modeCi.MODEJ.MODE)  +  d.LINK.km 

CALL  BINS_2(  i.MODE,  j.MODE,  f.MUF.MHzO.MODEJ.MODE),  - 
fMEIz  low,  fh5lz  hgh,  4, 4,  fMHz  modejbin,  • 

NFRQ) 

"  CALL  BINS_2(  i.MODE,  j_MODE,  a.MUF.dgCi_MODEJ.MODE).  - 

"  adgjow,  adg.hg^  4, 4,  adgjnode  bin,  N  ANG  ) 

"  CALL  BINS_2(  i.MODE,  j_MODE,  d.MOT.kmOlMODEJ.MODE,!),  - 

•  dkm.low,  dkm.hgh,  4, 4,  dan.mode.bm,  N.DST  ) 

CALL  BINS_3(  i.MODE,  j_MODE,  i.HOUR,  f.MUF.MHzO.MODEJ.MODE),  - 
fMHz  low,  fMHz  hgh,  4, 4, 24,  fMHz  mode  hr.bin,  - 
NFRQ) 

"  CALL  BINS.3(  i.MODE,  j.MODE,  i.HOUR,  a.MUF.dg(i.MODEJ_MODE),  - 

"  adgjow,  adg.hg^  4, 4, 24,  adgjmode.hr  bin,  - 

"  N.ANG) 

"  CALL  BINS_3(  i.MODE,  j.MODE,  i.HOUR,  d.MUF.km(i.MODEJ_MODE,  1),  - 

"  dkm  low,  dkm.hgh,  4, 4, 24,  dkm.mode.hr.bin, - 

-  N.DST) 


END  DO 
END  DO 

DOi.DISP-l,N.DISP 
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IF  (  f  LINK  MHz  GT.  f  DlSP.l)  - 

.AND.  f jLINK_MHz  .LE.  fJimit(LDISP,2) )  THEN 

fMHzJir  disp(i  HOUR,i  DISP)  -  {MHz^hr_disp(i_HOUR,i_DISP)  +  1.0 
iMHzIdisp(i_DiSP)  « fl^_disp(iJ)IS^  +  1.0 

END  IF 

END  DO 

"  Update  frequency-conditional  take-ofT  angle  mean  bins. 

LFRQ  -  INT(  (  f_HNK_MHz  -  fMHz Jow  )  /  d.fMHz  )  +  1 

IF  (  f_LINK_MHz  1£.  fMHzJow  )  THEN 

ANG  hr  frqO  HOUR,0.1)  -  ANG_hr  frq(i_HOUR,0,l)+ 1.0 
ANG.hr^frqO JIOUR,0,2)  *  ANG_hr_frq(i_HOUR,0,2)  +  a^^LINK^dg 

ELSE  IF  (  f.LINK_MHz  .GT.  fMHzJigh  )  THEN 

ANG.hr  firqCiJiOUR,N  FRQ+1,1)«=  ANG  hr  frq(i  HOUR,N  FRQ+1.1)- 

+  1.0 

ANGJirJirq(iJIOUR,N_FRQ+l,2)  -  ANGJir.firq(iJHOUR,NJFRQ+l,2)  - 
+  a.LINK.dg 


ELSE 

ANG  hr  frqO  HOUR,i  FRQ.l)  -  ANG.hr  frq(i_HOUR,i_FRQ,l)+  1.0 
ANG.hr Jrq(iJIOUR,i_FRQ,2)  =  ANGJir_frq(i_HOUR,i_FRQ,2)  +  a.LINK.dg 

END  IF 

RETURN 

END 


"  SUBROUTINE  k.VALUE 

N 

*  Purpose:  To  compute  the 'secant-corrected' effect  of  a  curved  ionosphere 


REAL*4  FUNCTION  k_VALUE(  D  ) 


A-61 


REAL*4  D,  D_8mp(17),  k;_smp(17),  m_k,  b_k 

DATA  D_smp/0.0,23 1 .9,3 18.8,376.8,492.8,608.7,782.6,985.5, 1246.4,- 
2800.0,3014.5,3205.8,3373.9,3536.2,3687.0,3884.1,- 
4000.0/ 

DATAk_sinp/1.0,1.0,1.0007, 1.0014,1.0035,1.0058.1.0111,1.0188,- 
1.0294,1.1047,1.1153,1.1259,1.1365,1.1471,1.1612,- 
1.1772,1.1871/ 

IF (D  ,LT. D  smp(l) ) THEN 
k,VALUE-k  smp(l) 

RETURN 

END  IF 

IF  (  D  .GE.  D_smp(17)  )  THEN 
k^VALUE  *  k_smp(17) 

RETURN 

END  IF 

DOk«2, 17 


IF  (  D_smp(k-1)  .LE.  D  .AND.  D  .LT.  D_smp(k)  )  THEN 

DUM_1  ■  D_smp(k)  -  Djnnp(k-1) 
mjc  «  (  k_smp(k)  -  k_si^-l)  )  /  DUMJ 

■  (  k_8mp(k-l)  *  D_smp(k)  -  k_smpi)  ♦  D_smp(k-1) )  /  DUM  1 
k_VALUE  -  mJc  ♦  D  +  b_k 
RETURN 

END  IF 

END  DO 

RETURN 

END 


"  SUBROUTINE  binFILO 

N 

"  Purpose:  Increment  density  function  bin  corresponding  to  input  value 


SUBROUTINE  BINS_0(  Vr,  Vr_low,  Vr_hgh,  bin__0,  N_bins  ) 
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REALM  bin_0(0:61) 

REAL*4  Vr,  Vr  low,  Vr  hgh,  d  Vr 
INTEGERMNjins 

LOGICAL*!  B_flag 

IF(Vr.LE.VrJow)THEN 
bin  0(0)*  bin  ()(0)-t- 1.0 
ELSE  ff  (  Vr  .GTrVr_hgh  )  THEN 

bin  0(N  bins+l)*bin  0(N  bins+l)+1.0 

ELSE 
ij)in  « 1 

B^flag*  JALSE. 

d_Vr  -  (  Vr^hgh  -  Vr  low  )  /  FLOAT(  N  bins  ) 

DO  WHILE  (T_bin  .LE”N_bins  .AND.  .NOT.  B_flag  ) 
Vr_min  *  Vr__Iow  +  ( ijbin  -  1  )  *  d_Vr 
Vrjnax  *  Vr_min  +  d  Vr 

IF  (  Vr  .GT.  Vr_inin  Vr  .LE.  Vrjnax  )  THEN 

bin_0(ij)in)  *  bin_0(i_bin)  +1.0 
Bjiag*.TRUE. 

ENDff 

i  bin  *  i  bin  +  1 
EJ®DO  ~ 

END  IF 

RETURN 

END 


SUBROUTINE  binFILl 

Purpose:  Increment  density  function  bin  corresponding  to  input  value 

SUBROUTINE  BINS_1(  i_l,  Vr,  Vrjow,  Vrjigh,  N_l,  bin_l,  N_bins  ) 

REAL*4bin_l(l;24,0:61) 

REAL*4  Vr,  Vr  low,  VrJigh,  d_Vr 
INTEGER*4  N^bins,  i_l,  N_1 

LOGICAL*!  B^flag 

IF(Vr.LE.Vr  low) THEN 

bin  l(i  l,03-bin  l(i  1,0) +1.0 
ELSE  (  Vr  .GT.  Vr_hgh )  THEN 

bin_l(i_l,N_bins+l)  *  bin_l(i_l,NJbms+l)  +  1.0 


FISF 
i_bin  ■  1 

B  flag  « .FALSE. 

"d  Vr  -  (  Vr  hgh  -  Vr  low  )  /  FLOAT(  N_bins  ) 
DOWHELE (T_bm  .LE”N_biiis  .AND.  .NOT.  BJlag ) 
Vr__iiiin  *  Vr^low  +  ( ijbin  -  1 )  *  d_Vr 
Vr  max»Vr  min-t-d  Vr 

IF  (  Vr  .GT.  Vr^min  .AND.  Vr  LE.  Vr^max  )  THEN 
bm_l(i_l,i_bin)  *  bin_l(i_l,i_bin)  +  1.0 
B  fiag«.‘n^UE. 

ENDff 

i  bin  » i  bin  +  1 
EMJDO  " 

END  IF 

RETURN 

END 


"  SUBROUTINE  binFIL2 

II 

”  Purpose:  Inaement  density  function  bin  corresponding  to  input  value 

SUBROUTINE  BINS_2(  i^l.  ij,  Vr.  Vrjow,  Vrjigh,  N^l,  N_2,  - 
binJZ,  NJbins  ) 

REAL*4  bin_2(l:4, 1:4,0:61) 

REALM  Vr,  Vrjow,  Vr_hgh,  d_Vr 
INTEGERM  NJbins,  i_l,  ijl,  N_l,  N_2 

LOGICAL*!  BJlag 

IF  (  Vr  .LE.  Vrjow  )  THEN 

bin  2(i  l,i  2,0)»bin_2(U,i_2,0)+ 1.0 
ELSE# (Vr.GT.  Vr  hgh)TH]^ 

bin  2(iJ,i_2,NJ)ins+l)  =  bin_2(i_l,i_2,N_bins+l)  +  1.0 

ELSE 
i  bin*  1 

B  flag*  JALSE. 

hyr  »  (  Vr_hgh  -  Vrjow  )  /  FLOAT(  NJ)ins  ) 

DO  WHILE  ( i  Jin  iE.  NJ)ins  .AND.  .NOT.  B_flag  ) 

Vr_niin  *  Vrjow  +  ( i J)in  -  1 )  *  d_Vr 
Vr_max  *  Vrjaain  +  d_Vr 

If7  Vr  .GT.  Vr_inin  .AND.  Vr  .LE.  Vr_max  )  THEN 

bin_2(i_l,r2,iJ>“)  “  bm_2(Ll.L2,iJ>in)  +  10 
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B  flag  -  .TRUE. 
ENDff 

i  bin  *  i_bin  +  I 
EM>DO” 

END  IF 

RETURN 

END 


SUBROUTINE  binFILS 

Purpose;  Increment  density  function  bin  corresponding  to  input  value 


SUBROUTINE  BINS_3(  i_l,  i  2,  i_3,  Vr.  Vr  low.  Vr_hgh.  N  1,  N  2,  - 
N_3,  bin_3.N_bins) 

REALM  binj(l:4, 1:4, 1:24,0:61) 

REAL*4  Vr,  Vr  low,  Vr  hgh,  d  Vr 
INTEGER*4  Njins,  LC ij,  O,  NJ2,  N_3 

LOGICAL*!  BJlag 

IF  (  Vr  .LE.  Vrjow  )  THEN 

bin  3(i  l,i  2,i  3,0) -bin  3(i_l,i  2,i_3,0)+1.0 
ELSE  ff  (  Vr  .GT.  Vr_ligh  )  THEN 

bin_3(i  l,i_2,i_3,N  bins+l)-bin  3(i_l,i_2,i_3,N  bins+l)+1.0 

ELSE 

i_bin-l 

B  flag -.FALSE. 

‘d_Vr  -  (  Vrjigh  -  Vrjow  )  /  FLOAT(  N  bins  ) 

DO  WHILE  ( ij)in  .LE.  N_bins  .AND.  .NOT.  B_flag  ) 

Vr_inin  -  Vrjow  +  ( i J)in  - 1 )  *  d_Vr 
Vr_max  -  Vr_inin  +  d_Vr 

IF  (  Vr  .GT.  Vr  min  .AND.  Vr  .LE.  Vr_max  )  THEN 
bin_3(iJ,r2,L3,Lbin)  -  bin_3(Ll,L2,L3,iJ)in)  +  1.0 
B_flag  =  .TTIUE. 

END  IF 

ij)in  -  i J)in  +  1 
END  DO 
END  IF 


RETURN 


END 


•  SUBROUTINE  PARMBN2 

m 

•  Purpose;  Eo^  parameter  values  from  UAG-23.NEW  formatted  files  into 
"  ^propriate  arrays  holding  all  hourly  values  for  a  given  month. 

SUBROUTINE  PARMBN2(  i_SITE.  LDAY,  i_HOUR,  i_CODE.  P_value,  Q_char,  - 
D_char  ) 

REAL*4  P_value 

INTEGER*4  i  SITE,  i_DAY.  i  HOUR,  i  CODE 
CHARACTER*!  CLchar.  D_char 

REAL*4  foF2,  fxF2,  fzF2,  M3000F2,  hF2,  hpf2,  foFl,  ficFl,- 
M3000F1,  hFl.  hF,  foE,  foE2,  hE,  hE2,  foEs,- 
fxEs,  fbEs,  fEs,  hEs,  foFldS.  iinin,  M3000Fld5,  hfldS,- 
fol,  fid,  finl,  12000,  n,  Ixxx 

"  CHARACTER*!  QfoF2,  Q£xF2.  Q£zF2,  QM3000F2,  QhF2,  QhpfZ,  QfoF!,  QfxF!,- 

-  QM3000F!,  QhF!,  QhF,  QfoE,  QfoE2,  QhE,  QhE2,  QfoEs,- 

"  QficEs,  QfbEs,  QfEs,  QhEs,  QfoF!d5,  Qfmin,  QM3000F!d5,. 

"  Qhf!d5,  Qfol,  Qfxl,  Qfinl,  QI2000,  QI,  Qbtxx 

"  CHARACTER*!  DfoF2,  I>fitF2,  Df2f2,  DM3000F2,  DhF2,  Dhpf2,  DfoF!,  DfxF!,- 
"  DM3000F!,  DhF!,  DhF,  DfoE,  DfoE2,  DhE,  DhE2,  DfoEs,- 

"  DfidEs,  DfbEs,  D£Es,  DhEs,  DfoF!d5,  Dfinin,  DM3000F!d5,- 

"  Dhfld5,DfoI,DfxI,DfinI,DI2000,DI,DIxxx 

H  0m*********0****»*0**m*******************m*******0*00)tt000**0*0*0titm 

”  Define  common  blodcs 

COMMON/IONSNP/  foF2(2,31,24),  M3000F2(2,31,24),  hF2(2,31,24),  - 
foF!(2,3!,24),  M3000F  1(2,3!, 24),  hF!(2,31,24),  - 
hF(2.31,24),foE(2,31,24),- 
hE(2.3!,24),  foEs(2,31,24),  fbEs(2,3!,24),  - 
£Es(2,3!,24),  hEs(2,3!,24) 

-  fiiiin(2,3!,24),fxl(2,3!,24) 

IF(i  CODEEQ.  !)THEN 
f<^i  SITE,!  DAY,i  HOUR)-P.va!ue 

•  (JfoF^i  SrrCi  DAY^i  HOUR)-(Lchar 

•  Df6F2Ci  SITE,!  DAY.i  HOUR)-D_char 
END  IF 


"  IF(LC0DE.EQ.2)THEN 

"  fxF2(i  SITE,i  DAY.i  HOUR)  »  P_value 

"  QfxF20_SITE,i_DAY.i,HO»JR)  *  Q_char 

"  DficF2(i  SITE,i  DAY.fHOUR)  =  D_char 

"  ENDIF 

"  IF(i  CODE.EQ.3)THEN 

"  f2F2(i  SnE,i  DAY.i  HOUR)-P  value 

"  Qf2F20_SITE,i_DAY,i  HOUR)-CLchar 

"  Df2F2(i  SITE,!  DAY,i_HOUR)  =  D  char 

"  ENDBF 

IF  ( i_CODE  .EQ.  4  )  THEN 
M3^F2(i_SITE,i  DAY.i  HOUR)  »  P_value 
"  QM3000F2(t  SITEii  DAY^i  HOUR)«Q_char 

"  DM3000F2(i  SITE,i  DAY.iJHOUR)  *  D  char 

ENDIF 

IF  ( i_CODE  .EQ.  5  )  THEN 
hF2(LSrrE.i  day,!  H0UR)-P  value 
“  QhF2(tSITE;i  DAVji  HOUR)  ^-CLchar 
"  DhF2(LSITE,i_DAY,i_HOUR)«D_char 
ENDIF 

"  IF  ( i_CODE  .EQ.  6  )  THEN 

"  hpK(i  SITE.i_DAY.i  HOUR)  »  P_value 

"  QhpF2(i  SITE,i_DAY,i  HOUR)  =  Q_char 

"  DhpF2(i  SrrE,i_DAY,i_HOUR)*D  char 

"  ENDIF 

IF(i  CODEEQ.7)THEN 
foFl(i  SITE,i_DAY,i_HOUR)«P  value 
"  QfoFKi  SITE,i_DAY,i_HOUR)»Q_char 

"  DfoFl(i_SITE,iJDAY,i_HOUR) »  D_char 

ENDff 

"  IF(i_CODEJEQ.8)THEN 

"  fxFl(i  SnE,i  DAY,iJHOUR)«P_value 

"  QfxFlO  SrrCi  DAY,i  HOUR)»Q_char 

"  DfxFl(LSITE,i_pAY,i_HOUR)«D_char 

"  ENDIF 

IF  ( i_CODE  £Q.  9  )  THEN 
M3000F10  SITE,!  DAY,i  HOUR)  «  Pjvalue 
•  QNDOOOFIO  SirCi  DA^i  HOUR)  =  Q_char 
DM3000Fl(i  SITE,!  DAY.i  HOUR)-D  char 


II 


END  IF 


IF(LCODE.EQ.  10)  THEN 
hFT(LSITE,i  DAY.i  HOUR)-?  value 
QhFl(i  SITE,i  DA^i  HOUR)-Q_char 
DhFl(i  SITE.rDAY,rHOUR)«D  char 
END  IF 

IF(i_CODE£Q.  11)  THEN 
hF(i  SITE,!  DAY,i  HOUR)-?  value 
QhF(LSIT^i  DAY^i  HOUR) -“Q_char 
DhF(i  SITE,!  DAY,rHOUR)-D  char 
END  IF 

IF  (  LCODE  .EQ.  12  )  THEN 
foE(LSITE,LDAY.i_HOUR)-?  value 
Qfc^LSITEii  DAY.i  HOUR) -l^char 
DfoE(i  SnE,rDAY,rHOUR)«D  char 
ENDBF 


IF  ( i_CODE  .EQ.  13  )  THEN 
foK(i  SITE,!  DAY,i  HOUR)-?  value 
QfoE2(i  SITeTi  DAYTi  HOUR)-CLchar 
DfoE2(i  SITE,!  DAY.i  HOUR)-D  char 
END  IF 

IF  (  LCODE  .EQ.  14  )  THEN 
h]^_SITE,i  DAY,i  HOUR)-?  value 
QhE(LSIT£;i_DA^i  HOUR) -"Q^char 
DhE(LSITE,i_DAY.rHOUR)  -  D^char 
ENDIF 


IF(i  CODE.EQ.  15)  THEN 
hE2(i  SITE,i_DAY,i_HOUR)»?  value 
QhE2a_SITE,i  DAY.i  HOUR)-CLciMr 
DhE2(i  SITE,iJDAY,rHOUR)  =  D_char 
ENDIF  " 

IF(i  CODE.EQ.  16)THEN 
foEs(i  SrTE,i_DAY,i_HOUR)  -  Fjvalue 
QfoEs(i_SITE,i  DAY,i_HOUR)  -  Q_char 
DfoEs(i  SITE.i  DAY.i  HOUR)-D_char 
ENDIF 


IF(i  CODE.EQ.  17)THEN 
fic^i  SnE,i  DAY,i  HOUR)  -  F.value 
QfxEs(i  SIT^i  DAY.i  HOUR)-Q_char 


“  DfxEs(i  SITE,!  DAY,i  HOUR)«D  char 

"  ENDIF 

IF(i  CODE.EQ.  18)  THEN 
£bB(i  SnE,i  DAY,i  HOUR)-P  value 
QfbEsO  SITE,!  DAYTi  HOUR) -"c^char 
DfbEsO  SITE,!  DAY,i  HOUR)-D  char 
ENDIF 

ff  (i  CODE.EQ.  19) THEN 
fEsO  SITE,i  DAY,i  HOUR)-P  value 
"  QfEsO.SrrEii  DA^i  HOUR)«Q_char 

"  DfEs(i  SITE,!_DAY,i  HOUR)-D  char 

ENDIF 

IF  ( i_CODE  .EQ.  20  )  THEN 
hEs(i_SITE.i  DAY.i  HOUR)  -  P_value 
"  QhEs(i_SITE;i  DAY^i  HOUR)  -  Q_char 
"  DhEs(i  SITE,!  DAY.i  HOUR)»D  char 
ENDIF 


IF(i  CODE  .EQ.  21)  THEN 
foFldSO  SrnE,i  DAY,i  HOUR)  =  P  value 
QfoFld5(i  SniE,i  DAY.i  HOUR)-'^char 
DfoFld5(i  SITE,i  DAY.i  HOUR)-D  char 
ENDIF 

IF  (i  CODE  £Q.  22  )  THEN 
fiiiin(i_SITE,i  DAY,i  HOUR)  -  P_value 
Qfinin(LSIT^i  DAYJ  HOUR)  =  Q_char 
Dfinin(i  SITE,i_DAY,rHOUR)  »  D  char 
ENDIF 

IF  ( i_CODE  £Q.  23  )  THEN 
M3000Fld5(i_SITE,i_DAY,iJIOUR)  »  Pjvalue 
QM3000Fld5(i_SITE,i_DAY,i_HOUR)  '"^char 
DM3000Fld5(i_SITE,i_DAY,i_HOUR)  =  D_^char 
ENDIF 

IF(i  CODE  .EQ.  24)  THEN 
hFld5(i  SITE.i  DAY,i_HOUR)  «  P_value 
QhFldSO  SIT^J>AY,i_HOUR)-’CLchar 
DhFld5(rsITE,i_DAY,i_HOUR)  =  D_char 


ENDIF 

ff  (i  CODE  .EQ.  25) THEN 
foI0_SITE,LDAY,i_HOUR)  -  P.value 


r 


QfoI(i  SITE,i  DAY,i  HOUR)  =  CLchar 
Dfol(i  SITE.i  DAY,i  HOUR)  =  D_char 
END  IF 

IF(i  CODE  .EQ.  26)  THEN 
fxIO  SITE,iJDAY,i  HOUR)«P  value 
QfidO  SITE,!  DAY,!  HOUR)-^char 
D£xI(rsiTE,rDAY,i  HOUR)-D  char 
ENDF 

IF  ( i  CODE  .EQ.  27  )  THEN 
fini(LSnE,LDAY,i  HOUR)-?  value 
Qfml(i  SITE,!  DAY.i  HOUR)«‘^char 
DfmlO  SITE,!  DAY,rHOUR)-D  char 
ENDIF 

IF  ( i_CODE  EQ.  28  )  THEN 
I2000(i  SITE,!  DAY,i  HOUR)-?  value 
012000(1  SITE,!  DAyTi  HOUR)-‘^char 
DI2000(tSITE,i_DAY.i_HOUR)  -  D_char 
ENDIF 

IF(i  CODE  .EQ  .  29)  THEN 
no  SITE,!  DAY,i  HOUR)-?  value 
QI(i  SITE,!  DAY,i  HOUR)- ^char 
DIO_SITE,i_DAY,i  HOUR)-D  char 
END§ 

IF  ( i_CODE  EQ.  30  )  THEN 
IxxxO  SITE,i_DAY,i  HOUR)-?  value 
QLcxxO  SITE,!  DAY,!_HOUR) -TLchar 
DIxxxO  SITE,!  DAY,!  HOUR)-D_char 
ENDIF 

RETURN 

END 


*  Functicm  ROOT 

M 

*  Purpose:  the  value  of  Qnun  given  the  functions 

*  ?iO  and  dPrOdt  and  Function  PrMAX,  vduch  fmds  the 
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maximum  value  of  PrO  given  Q  and  the  time  at  which 
this  maximum  occurs,  i.e.,  t_MAX 


REAL*4  FUNCTION  ROOT(  a_LOW_rd,  a_HGH_rd,  r_km.  As_rd.  INTV,  ITRN ) 

REAL*4  theta^diff,  r_km(2) 

LOGICAL*!  INTV,  ITRN. D_^flag 

PARAMETER(  ITMAX  »  100,  EPS  «  3.0e-8  ) 

ROOT  ERR  =  0.00001 
INTV -.TRUE. 

ITRN  -  .TRUE. 

D_flag-. FALSE. 

A-a_LOW_rd 

FA  -  theta_difi(  A,  aJLOW_rd,  a_HGH_rd,  rjcm,  As_rd ) 

B  =  a_HGH_rd 

FB  =  Aeta_diffl(  B,  aJLOW_rd.  aJIGH^rd,  r_km,  As_rd ) 

IF(FB*FA.GT.  0.0)THEN 
INTV -.FALSE. 

IF  (  FA  .LT.  0.0  .AND.  FB  XT.  0.0  )  THEN 
ROOT«a_LOW  rd 
RETURN 
END  IF 

IF  (  FA  .GT.  0.0  .AND.  FB  .GT.  0.0  )  ROOT  -  a_HGH_rd 
RETURN 
END  IF 
FC-FB 

DO  ITER- 1,  ITMAX 

IF(FB*FC  .GT.  0.0)'IHEN 

C-A 

FC-FA 

D-B-A 
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E-D 


END  IF 

IF  (  ABS(  FC  )  XT.  ABS(  FB  )  )  THEN 

A*B 

B-C 

C-A 

FA»FB 

FB-FC 

FC-FA 

END  IF 

TOLl  =  ::.0  •EPS  *  ABS( B )  +  0.5  •  ROOTJERR 
XM«0.5*(C-B) 

IF  (  ABS(  XM  )  XE.  TOLl  .OR  FB  XQ.  0.0  )  THEN 

IF(B.LT.a  LOW  rd  +  ROOT  ERR)THEN 
ROOT  «  a  LOW  rd 
ELSE 

ROOT  =  B 
END  IF 

RETURN 

ENDff 

IF  (  ABS(  E  )  .GE.  TOLl  J\ND.  ABS(  FA  )  .GT.  ABS(  FB  )  )  THEN 

S  =  FB/FA 

IF(AXQ.C)THEN 
P  =  2.0*XM*S 
Q-l.O-S 
ELSE 

Q-FA/FC 

R-FB/FC 

P-S*(2.0*XM*Q*(Q-R)- 

.(B-A)*(R.1.0)) 

Q»(Q-1.0)*(R-1.0)*(S-1.0) 

END  IF 
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IF(P  .GT.  0.0)Q--Q 
P-ABS(P) 

Y-MIN(3.0*XM*Q-ABS(TOL1  •Q).  ABS(E*Q)) 
IF(2.0*P.LT.  Y)THEN 
E-D 
D-P/Q 
ELSE 
D«XM 
E-D 
END  IF 

ELSE 

D-XM 

E-D 

END  IF 

A-B 

FA-FB 

IF  (  ABS(  D  )  .GT.  TOLl  )  THEN 
B-B  +  D 
ELSE 

B-B  +  SIGN(TOLl,XM) 

END  IF 

FB  -  theta_difil[  B,  a_LOW_rd,  a_HGH_rd,  rjcm,  As_rd  ) 

END  DO 

ITRN-JFALSE. 

ROOT-B 

RETURN 

END 


REALM  FUNCTION  thcta_diff(  a_rd,  a_LOW_rd,  a_HGH_rd,  r^km,  As_rd  ) 
REALM  DUM,  r_km(2),  As_rd,  Re_km,  a_km,  theta_rd(2),  alfa_rd(2) 
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REALM  Hearth  km 


INTEGERMLHOP 

COMMON  /MISC/  Re_km.  Rearth.km 

alfa_rd(])  •  a_rd 
alfa.rd(2)  ■  As_rd  -  a_rd 

DOi_HOP-l,2 

DUM-rJan(i  HOP)**2  +  Rearth  km**2 

DUM  -  DUM  -  2.0  •  r  km(i  HOP)  •  Rearth  km  •  COS(  alfa_rd(i  HOP)  ) 
a_km-SQRT(DUM) 

DUM  -  r,km(i  HOP)**2  -  a  km**2  -  Rearth_km**2 

DUM-DUM/(2.0*Rearth.km*a  km) 

IF  (  DUM  .GT.  1.0  )  DUM  -  1.0 
IF  (  DUM  .LT.  -1.0  )  DUM  *  -1 .0 
theta_rd(i_HOP)  -  ASIN(  DUM  ) 

END  DO 

theta^diff  « thetajrd(l)  -  theta_rd(2) 

RETURN 

END 
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APPENDIX  B 


MEASUREMENT-PREDICTION  COMPARISON  OF  EMPIRICAL 
FREQUENCY  AVAILABILITY  VALUES 


This  appendix  contains  comparison  of  empiricaliy-derived  frequency  availability  on  trans-polar 
two-hop  HF  radio  links  with  ASAPS  and  lONCAP  predictions.  The  empirical  data  was  derived  from 
UAG-23-formatted  ionosonde  data  measured  at  Dikson,  Cape  Schmidt,  and  Loparskaya  in  Russia  and 
Resolute  Bay  in  Canada  in  1991.  The  1991  year  corresponded  to  mid-range  sunspot  numbers  from  133 
to  172.  The  following  table  provides  a  description  of  the  plots  as  well  as  an  index  to  locate  specific 
plots. 

TABLE  B.1  INDEX  TO  FREQUENCY  AVAILABILITY  PLOTS 


FOT 

MUF 

HPF  1 

LINK 

RANGE 

MONTH 

PAGE 

MONTH 

PAGE 

MONTH 

PAGE 

JAN 

B-3 

JAN 

B-4 

JAN 

B-5 

DIKSON  TO 

FEB 

B-6 

FEB 

B-7 

FEB 

B-8 

RESOLUTE 

7100  KM 

MAR 

B-9 

MAR 

MAR 

B-11 

BAY 

APR 

B-12 

APR 

APR 

B-14 

MAY 

B-15 

MAY 

MAY 

B-1 7 

JUN 

B-18 

JUN 

B-19 

JUN 

B-20 

CAPE 

MAR 

MAR 

B-23 

MAR 

B-24 

SCHMIDT 

9100  KM 

JUL 

JUL 

B-26 

JUL 

B-27 

TO 

SEP 

B-28 

SEP 

B-29 

SEP 

B-30 

LOPARSKAYA 

DEC 

B-31 

DEC 

B-32 

DEC 

B-33 

B-1 


TABLE  B.2  INDEX  TO  FREQUENCY  AVAHABILITy  PLOTS 


FOT 

MUF 

HPF  1 

1  LINK 

RANGE 

MONTH 

PAGE 

MONTH 

PAGE 

MONTH 

PAGE  1 

^  ^  J;  ^ 

y<^».yyy  ’’W”TO 

•6  ' 
o\a'^ 

'  N  a' 

^  V.  ^  V  «>'  S/ 

>  '  ,  'i?  - 

>V  S  '  A 

~  yx\  ' 

„-A:  .rM:r  '  '.'  ib?  .  v'- 

^  ••  aNv  ••  ■^<’'  <<  $  ‘■s>-  *».  4 

CAPE 

SCHMIDT 

TO 

LOPARSKAYA 

9100  KM 

MAR 

JUL 

SEP 
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B-22 

B-25 

B.28 

B.31 

MAR 

JUL 

SEP 

DEC 

B-23 

B-26 

B-29 

B-32 

MAR 

JUL 

SEP 

DEC 

B-24 

B-27 

B-30 

B-33 

B-2 


FOT  COMPARISON  FEB  1991  SSN:  133  (AS APS) 
RANGE:  7100  KM  DIKSON  &  RESOLUTE  BAY 


lONOSONDES  — lONCAP  ASAPS 


FOT  COMPARISON  MAR  1991  SSN:  181  (ASAPS) 
RANGE:  7100  KM  DIKSON  &  RESOLUTE  BAY 
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SdVSV  dVONOI  SaONOSONOI 


MUF  COMPARISON  MAR  1991  SSN:  181  (ASAPS) 
RANGE:  7100  KM  DIKSON  &  RESOLUTE  BAY 
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lONOSONDES  —  lONCAP  ASAPS 


HPF  COMPARISON  MAR  1991  SSN:  181  (ASAPS) 
RANGE:  7100  KM  DIKSON  &  RESOLUTE  BAY 


B-11 


lONOSONDES  — lONCAP  — ASAPS 


FOT  COMPARISON  APR  1991  SSN:  172(ASAPS) 
RANGE:  7100  KM  DIKSON  &  RESOLUTE  BAY 


SdVSV  dVONOI  —  SaONOSONOl 


MUF  COMPARISON  APR  1991  SSN:  172(ASAPS) 
RANGE:  7100  KM  DIKSON  &  RESOLUTE  BAY 
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lONOSONDES  — ^  lONCAP  ASAPS 


HPF  COMPARISON  APR  1991  SSN:  172(ASAPS) 
RANGE:  7100  KM  DIKSON  &  RESOLUTE  BAY 
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lONOSONDES  •—  lONCAP  ASAPS 


FOT  COMPARISON  MAY  1991  SSN:  143(ASAPS) 
RANGE:  7100  KM  DIKSON  &  RESOLUTE  BAY 


CM 


B-15 


lONOSONDES  — ^  lONCAP  ASAPS 


MUF  COMPARISON  MAY  1991  SSN:  143(ASAPS) 
RANGE:  7100  KM  DIKSON  &  RESOLUTE  BAY 


CM 


B-16 


lONOSONDES  — lONCAP  ASAPS 


FOT  COMPARISON  JUN  1991  SSN:  133  (ASAPS) 
RANGE:  7100  KM  DIKSON  &  RESOLUTE  BAY 


B-18 


lONOSONDES  — ^  lONCAP  — ^  ASAPS 


MUF  COMPARISON  JUN  1991  SSN:  133(ASAPS) 
RANGE:  7100  KM  DIKSON  &  RESOLUTE  BAY 
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lONOSONDES  — ^  lONCAP  ASAPS 


TABLE  INDEX  TO  FREQUENCY  AVAILABILITY  PLOTS 


FOT 

MUF 

HPF  1 

LINK 

RANGE 

MONTH 

PAGE 

MONTH 

PAGE 

MONTH 

PAGE 

JAN 

B-3 
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FEB 
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HPF  COMPARISON  MAR  1991  SSN:  181  (ASAPS) 
RANGE:  9100  KM  CAPE  SCHMIDT  &  LOPARSKAYA 


B<24 


lONOSONDES  lONCAP  — ^  ASAPS 


FOT  COMPARISON  JUL1991  SSN:  150  (ASAPS) 
RANGE:  9100  KM  CAPE  SCHMIDT  &  LOPARSKAYA 
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lONOSONDES  lONCAP  — ASAPS 


HPF  COMPARISON  JUL1991  SSN:  150(ASAPS) 
RANGE:  9100  KM  CAPE  SCHMIDT  &  LOPARSKAYA 
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FOT  COMPARISON  SEP  1991  SSN:  148  (AS APS) 
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